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t3.  ABSTRACT 

.  in^xiinwdity  (ra  ^  1)  inventory  control  model  with  periodic  review  is  fonmia- 
t^  ^  studied  vAv^  a  Sredic  stationary  ordering  (d ,  S)  policy  in  in'  is  usedi  It 
is  ^SOTed  that:  (1)  the  ordering  decision  at  the  b^ihiring  of  eacdi  of  a  sef^ehce 
Of  p^iods  of  time  is  affected  ly  a  single  set  up  cost,  a  linear  variaJole  ordering  : 
cost  and  holding  and  shortage  costs;  (2)  tbe  demands  for  the  items  in  each  period 
.identically  and  independently  distributed  continuous  randem  variables;  (3)  the 
delivery  of  orders  is  immediate;  and  (4)  ccmp.lete  backlogging  of  infilled  demands 
is  ^iCTved.  Next  the  model  is  analyzed  by  minimizing  the  expression  for  the 
stationary  total  e^^obed  cost  per  period.  Tne  set  of  simultanecws  equations  used 
to  determine  the  optimal;  policy  parameters  are  restated  in  terms  of  a  real  valued 
function  in  iP.  Finally,  the  optimization  analysis  is  restricted  to  the  special 
cd£K2  of  a  ta^xj-ccnrnoditjT  invent»^rv  control  problem  idiere  tlie  demands  for  tlie  items 
obey  the  exponential  distriJ ,  and  the  holding  and  shortage  costs  are  linear. ( 
^:c  integral  cquatioji  v;hich  ii.  used  in  eolvii'jg  for  tliu  optical  policy  parameters 
is  converted  into  a  linear  hifpcrtolic  parti.al  differential  equation  of  the  second 
order  with  Loundary  condit.ions.  Ihis  boundury  value  problem  is  solved  analytically. 
The  so3utio:i  is  then  used  to  dotexirdne  the  cj^jLinal  polici’  parameters.  In  a 
;  numer-ical  cxan?)lc  a  computer  program  is  deve].o£32d  to  determine  the  optimal  policy 
parameters. 
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ABSTRACT 


« 

An  m-ccnrwdity  (m  ^  1)  iir/entory  cxAtrol  model  with  periodic  review 
is  fomulated  and  studied  when  a  dyadic  stationary  ordering  (a,  S)  policy 
in  is  used.  It  is  assuned  that:  (1)  the  ordering  decision  at  the 
beginning  of  each  of  a  sequence  of  periods  of  time  is  affected  a  single 
set  cost,  a  linear  variable  ordering  cost  and  holding  and  shortage  costs; 

(2)  the  demands  for  the  items  in  each  period  are  identically  and  independ^tly 
distributed  continuous  randan  variables;  (3)  the  delivery  of  orders  is 
imnediate;  and  (4)  coiplete  backlogging  of  unfilled  demands  is  allowed. 

Next  the  model  is  aneilyzed  by  minimizing  the  esqsression  for  the  stationary 
total  ejtpected  cost  per  period.  The  set  of  simultaneous  equations  used  to 
determine  the  optimal  policy  parameters  are  restated  in  terms  of  a  real 
valued  function  in  eP*.  Finally,  the  c^timization  analysis  is  restricted 
to  the  special  case  of  a  two-comodity  inventory  control  problem  where  the 
demands  for  the  items  cb^  the  e^qxmential  distribution,  and  the  holding 
and  Portage  costs  are  linear.  The  integral  equation  which  is  used  in 
solving  for  the  optimal  policy  parameters  is  converted  into  a  linear  hyper¬ 
bolic  partial  differential  equation  of  the  second  order  with  boundary  con¬ 
ditions.  This  boundary  value  problem  is  solved  analyticaTly.  The  solution 
is  then  used  to  determine  the  optimal  policy  parameters.  In  a  numerical 
exanple  a  ccmputer  program  is  develc^jed  to  determine  the  optimal  policy 
parameters. 
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CHAPTEK  1 


INERDDUCnON 

1«1  Introduction 

Ohe  replenishment  and  cxintxol  of  inventories  is  a  cannon 
practice  in  the  fields  of  business,  eoononics  and  nanagement  in  general. 
Sijiply,  an  inventory  is  a  stock  of  physical  goods  which  are  held  or 
stored  for  future  use,  sale  or  production.  Inventory  prdolems  are 
integrated  in  nature.  involve  production,  schedulirjg,  distri¬ 

bution  of  oanmodities,  or  oorbination  of  all.  In  the  past  twenty  years, 
ituch  work  has  been  done  to  fcoulate  and  a^ply  the  proper  nathanatical 
tools  for  determining  economic  decision  rules  in  managing  the  various 
inventory  ^sterns.  Viewed  as  decision  problems,  the  analysis  of  sudi 
systons  essentially  proceeds  in  three  steps:  (1)  the  formulation  of  a 
'  model  e^^ressing  a  set  of  r(!lations  between  a  set  of  decision  variables 
\dx>se  values  are  to  be  determined;  (2)  the  establishment  of  an  objective 
function  of  the  decision  variables  so  as  to  minimize  or  maximize 
an  eoonanic  objective;  and  (3)  a  numerical  method  to  determine  the  values 
of  the  decision  variably  that  optimize  the  objective  function. 

In  this  introduction  we  shall  restrict  ourselves  to  models  in 
which  decisions  are  made  at  the  beginning  of  each  of  a  nuriber  of  equal 
periods  of  time,  and  in  vhich  the  demands  are  independently  and  identi-' 
cedly  distributed  random  variables  from  period  to  period.  Ihe  decision 
at  the  beginning  of  each  period  is  affected  by  a  fixed  set  vp  cost,  a 
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variable  order JLig  cost  and  a  holding  and  shortage  costs.  Backlogging 
and  imnediate  delivezy  of  orders  are  alleged. 

TVo  principal  distinct  approaches  have  been -used  in  analyzing 
such  inventory  models  both  in  theory  and  practice.  In  the  first  e^^jroach, 
the  system  is  viaited  as  a  multistage  decision  process  and  the  technique 
of  dynamic  prograimdng  is  enployed  in  finding  the  cptinal  policy  that 
minimizes  the  total  ejpected  cost  over  the  duration  of  the  process,  lb 
ensure  that  the  optimal  policy  is  of  the  "sinple  type"  certain  restric¬ 
tions  are  imposed  on  the  cost  factors  that  affect  the  decision  rules. 

More  precisely/  vhen  the  ordering  cost  is  linear  and  when  the  sum  of 
holding  oiid  shortage  costs  L(x)  is  convex/  the  optimal  policy  is  of  the 
"simple  type."  Hewever/  vhen  the  duration  of  the  process  is  infinite, 
a  secxjnd  approach  is  often  used:  A  "simple  type"  ordering  policy  is 
chosen  to  be  used  in  each  period.  Under  the  chosen  simple  policy  the 
sequence  of  inventory  levels  at  the  beginning  of  each  period  fonns  a 
Markov  process/  vhose  statiextary  behavior  can  be  analyzed.  Here/  instead 
of  optimizing  the  total  cost  over  the  horizon  period,  vhich  is  meaningless/ 
we  minimize  the  stationary  total  e:pecrf;ed  cost  per  period.  Ihe  optimal 
values  of  the  decision  variables  are  called  the  steadi-state  or  sta¬ 
tionary  solutiv-^ns  to  the  inventory  problem.  Seme  of  the  advantages  of 
the  stationary  approach  ever  the  dynamic  programming  approach  are: 

1.  Ihe  stationary  c^proach  provides  us  with  information 
daout  the  dependency  of  the  cptimal  policies  on  the  many 
parameters  involved  in  the  model  and  about  sensitivity  of 
costs  as  function  of  the  policies,  while  this  information  is 
not  provided  by  following  the  ^amic  programming  approach. 
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2»  As  the  nurbcr  of  oonmoditiGS  increases,  the  nurrber  of  the 
state  variables  increases  drastically,' and  the  use  of 

I 

dynamic  progranming  techniqjes  in  obtaining  optijnal  policies 
for  the  inventory  problem  will  be  time-consming  even  v^en 
using  the  fastest  ootputers.  ^ 

!?requently,  for  Icw-oost  high  turriiover  items,  optimal 

i 

policies  are  not  really  required.  It  is  sufficient  to  v 
obtain  simple  analytical  approximations  to  optimal  policies , 
and  these  are  obtained  fran  the  stationary  ^roach. 

Bie  main  ooncem  of  this  research  v;ill  be  to  find  the  stationary 
solutions  to  the  m-oonraodity  (m  >  1)  inventory  problem.  Per  the  one 
oomtDdity,  the  theory  of  inventory  has  studied  the  stationary  behavior 
of  the  ^stem  and  provided  optimal  decision  rules  for  inventory  rep?.en- 
ishment  binder  a  linear  ordering  cost  and  a  convex  holding  and  shortage 
costs  function.  Ihe  results  of  the  caie'^ocnnDdity  inventory  system  are 

^  t 

valid  for  multi-oonmoclity  inventory  ^tem  v^ere  each  contiodity  is 
treated  independently. 

:  I 

j 

However,  in  practice  this  independency  assunption  is  unoermon 
and  it  is  more  realistic  to  incorporate  the  interdependency  in  the  set 
tp  order  cost  of  the  oormodities  vhich  arises  in  the  nature  of  the 
follcwing  types  of  problems. 

'  i 

1.  Ordering  different  comodities  at  the  same  time  frem  a 
octimon  vendor  may  reduce  the  procurement  ordering  costs, 
thus  incurring  one  ordering  cost. 

2.  .Several  warehouses  used  as  storage  and  distribution  points 
are  supplied  by  a  single  factory.  Whenever  an  order  to 
rcploiish  a  given  oernnodity  stored  by  the  warehouses  is 
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made  to  the  factory/  a  fixed  set  tp  cost  is  incurred 
v^ethcr  the  warehouses  order  individually  or  simultaneously, 
dhe  configuration  of  this  system  is  shewn  in  Fig.  (1.1). 

Bie  fewo.prdalems  are  similar  in  structure  and  it  is  practical 
to  t^e  advantage  of  the  reduced  procurement  cost  in  controlling  the 
inventory  ,9^  the  oemnodities.  Furthermore/  if  one  can  solve  the  multi- 
ocranodity  problem/  tlien  the  multi-warehouse  problsn/  as  given  previously/ 
is  solved. 

Hie  procurement  cost  of  ordering  £  ^  ;0  units  is  assumed  to  be 

I 

equal  to  a  set  cost  plus  a  linear  purchase  cost 


C  2  ®  (^1 1  ^2  9***9  ®  ) 

“*  m 


m 

vAiere  ^  i  =  1/.../  m/  is  the  unit  purchasing  cost  v''f  item  i. 

Uie  fixed  set  rp  cost  is  incurred  only  if  an  order  is  made.  Thus  if 

I 

£  0  units  are  ordered/  the  fixed  cost  will  be  a,  function  of  z^,  say  K(£) 

In  genial  K(£)  will  take  as  many  different  values  as  the  notiber  of 

] 

I  alternative  different  ordering  decisions.  For  the  two  conmodity  problem 


K(2)  = 


Kl 

if  Zj  >  0  and  Z2 

K2 

if  Zj  «  0  and  Z2 

K 

if  Zx  >  0  and  Z2 

0 

if  Zi  =  0  and  Z2 

vdiere  Kj/  K2  ahd  K  are  all  nonnegative/  and  the  inequality  max  (Ki/  K2) 
^  K  ^  Kj  +  K2  is  satisfied.  ThrtS/  the  procurement  cost  is  given  ly 

:  *P  '  ! 

K(2)  +  C  2 


ii 

■-A 
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It 2  Literature  Review 

We  shall  first  sxirvey  the  literature  on  the  steady-state  or 
stationary  solutions  to  the  one-cxxnnnDdity  inventory  systan  operating 
under  a  (s,S)  policy;  next  we  shall  elaborate  on  the  work  done  on  the 
stationary  solutions  to  the  m-oorarodity  inventory  system  operating 
under  a  (o,S)  policy. 

Ibr  the  one-conmodity  problem  we  shall  consider  a  dynamic 
inventory  model.  At  the  beginning  of  each  pariod,  a  decision  is  made 
to  order  up  to  S  if  the  stock  falls  below  s,  otherwise,  nothing  is 
ordered.  The  demands,  in  successive  periods  (i  =  1,  2,  .  .  .) 

are  assumed  to  be  described  by  a  sequence  of  non-negative  random 
varidiles,  independently  and  identically  distributed  with  a  joint 
density  function  Delivery  of  orders  is  iimtediate  and  total 

backlogging  of  unfilled  donands  is  assumed.  The  ordering  decision  in 
each  period  is  affected  by  a  set  i^)  cost  K,  a  linear  purchase  cost  cz 
vhere  z  is  the  quantity  ordeired  and  c  the  unit  cost,  a  holding  cost  for 
carrying  inventory,  and  a  shortage  cost  for  not  meeting  d(2mands.  L(x) 
represents  the  expected  inventory  holding  and  shortage  costs  for  being 
in  stock  X  at  the  beginning  of  the  period. 

Ibllowing  the  pioneering  work  of  Arrow,  Harris,  Marchak  fl], 
Karlin  [2]  studies  the  steady  state  for  this  nodel  operating  under  a 
(s,S)  policy.  With  this  ordering  policy,  the  sequence  of  inventory 
levels  Xi,  X2,...,  vhere  X^^  is  the  inventory  level  at  the 

end  of  the  ith  period,  forms  a  Markov  proc:ess.  Karlin  shows  that 
the  distribution  of  coverages  in  the  sense  of  distributions  to  a 
stationary  distribution  with  the  density  f  (•)  given  by 
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\l>(s  -  X) 
i  +  y(s  -  s) 


for  s  <  X  <  S 


f(x) 


/S-s 

♦  (s  -  s)  +  I  .  i;/(t)  ^{S  -  X  -  t)  dt 

- ryyjs  l  b) -  for-<x<  s 


(1.1) 


v^eie 

♦(X)  «  I  (X) 

ml 

and 

¥(x)  =  I  (x) 

n=l 

♦  (•)  and  {•)  are  the  vrell-kncwn  syirbolic  notations  for  the 

fold  convolution  of  ({»(•)  and  $(•)»  respectively ►  Ihe  stationary 
distribution  is  clearly  independent  of  ai^  e(xncmic  consideration. 
By  iji|)osing  the  assuned  cost  stxuc:ture  on  the  process  the  total 
e}qpected  cost  per  period  denoted  by  g(s,  S-s)  is  given  as 


K  +  L(S)  +  L(S  -  X)  t|;(x)  dx 
g(8,  S  -  8}  - - - +  CM 


(1.2) 


vihexe  Q  «  S  -  s  and  m  ®  E(D) »  E(»)  is  the  expectation  operator. 


;  ’i  -  - 
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Greenberg  [5]  obtains  fjj(*)#  the  transient  soluticn  for  the 
(s,S)  inventory  system.  Using  the  generating  function  approach, 
the  stationary  distribution  f  (•)  is  given  by  the  standard  Abelian 
theorefn 


lim  ••  _ 

f  (x)  =  (1  -  z)  I  z*^  (x) 

z^l  n=l  ^ 


vhich  results  in  (1.1) .  lglehai±  [7]  shows  that  the  optimal  values 
of  s*  and  Q*  minimizing  g(s,  S  -  s)  are  given  as  the  solution  to  the 
following  equations 


L'(Q*.+  B*)  +  L'(Q*  +  -  x)  *(x)  dx  »  0 

K  +  L{Q*  +  s*)  +  j  L(Q*  +  s*  ~  x)  ^(x)  dx 
L(s*)  - - 2 - - - 


1  +  ♦(Q’) 


'vhere  L'(«'  is  the  derivative  of  L(*). 

Another  approach  in  studying  the  (s,S)  inventory  model  is  followed 
by  i>;vazlian  [11].  He  defines  to  be  the  total  expected  cost 

of  operating  the  intern  for  n  periods  vhen  no  order  is  placed 
initially,  i.e.,  s  5  x  <  S,  h^(x)  to  be  the  total  esqpectad  cost  of 
operating  the  ^tem  for  n  periods  vhen  an  order  is  placed  initially 
i.e.,  -•  <  X  <  s.  Then  fjj(x)  and  h^(x)  satisJ^  the  foUowing 
functional  equations 


I 


rx-s 


fjj(x)  =  L(x).  +  I 


+  j*  -  t)  ((i(t)  dt  ,  s  £  X  <  S  (n  ^  1) 


x-s 


,  /S-s 

h^(x)  =  K  +  c(g  -  x)  +  L(S)  + 


+  r  1<S  -  t)  ♦(t)  dt 
Js-s  ' 


-«  <  X  <  s  (n  >  1) 


Starting  wi*th  these  basic  functjional  equations  and  using  Howard's 

i  ‘ 

approadi  [6]  ^e  deduces  the  stationary  total  ejqoected  cost  p^iod 
as  given  by  (1.2) .  Ihe  optiinal  values  foif  s*  and  Q*  «  S*  -  s*  :that 

ndnimize  the  total  ejspected  cpst  per  period  are  given  as'  the  solvftion^ 

'  •  .  !  .1 

to  the  set  of  equations  ' 


M(s*,Q*)  =  K  and 


aM(s*yX) 

8x 


«  0 


xK)* 


lithcre 


Um(s*,x))  ==  Hl(s*)  ~  L(s*  •  xn 


the  operatpr  L  {  }  is  the  Laplace  transform  operator. 


£  *' 

P  ( 


S0/eral  cxarputationcil  methods  for  obtaining  values  of  s*  and  Q*  ate 
given  and  elaborated  later  by,  Sivazlian.  ^  for  the  case  of  gamna 
distributed  denuuid  and  linear  holding  and  shortage  costs,  a  diiicn*- 

sibnal  analysis  is  carried  to  reduce  the  nvnber  of  v^iables.  A 

•  '  '  '  .  ’ 

numerical  inversion  technique  of  the  Laplace  transform  e3q}ressions 

^  •» 

u$^  Gaussian  quadrature  is  used  to  solve  for  s*  and  Q*.  For  large 

•  ,  *  I 

values  of  Q,  Sivazlian  [11]  and  Roberts  [3]  shext^  t^t  s*  is  the 
solution  to  the  equation 

I  ‘  ■  ■  !  .  . 

'  f.  Cl'-  ♦Ml  d»  >  (I  (1+ 23  + 

>  <  » 

'  *  - _  1  .  ' 

and  Q  =  ’2K^'/h  '  variance  of 

the  randem  variable  D,  h  is  the  unit  holding  cost,  and  p  the  unit 
shortage  cost*  Sivazlian  [11]  also  shms  that  for  small  Q,s*  and 
Q*  ar?  the  solutions  for  the  follcx^g  appreodmate  system  of 


equations 


L'(s*)  >:  -  p 


L«(8*)  « 


*.  iK  .  2K 


0*7 -^5* 


'it  is  iaportant  to  note  that  an  integral  equation  of  the 
reneval  t^pe  in  one  dimension  is  solved  by  each  of  the  authors  in 

I  ♦  X 

carrying  their  ana?.^is;  In  Scarf's  work  [2] ,  the  density  function 

*  ,  I 

f  (•')  as  given  by  (1.1)  is  a  solution  bf  a  renewal  type  integral 


^See  .additional  referenoes> 


11 


equation.  Sivazlian  [11]  in  finding  the  total  expected  cost  per 
period  solves  a  renewal  type  integral  equation. 

For  the  case  v^en  the  demand  is  discretely  distributed, 
Veinott  and  Wagner  [16]  develop  the  nvethod  of  stationeiry  or  renewal 
analysis  for  the  oonputation  of  an  optimal  (s,S)  policy.  The  total 
e3q)ected  cost  per  period  is  givai  by 


g(s,  S  -  s)  = 


S-s 

K  +  L(S)  +  X  L(S  -  j)  t^j) 

j=l 

1  +  '}'{S  -  s) 


+  cy 


%Aiich  is  a  discrete  version  of  (1.2) .  Ihe  algorithm  provided  for 
searching  the  (s,S)  policies  to  find  one  that  minimizes  g{s,S  -  s) 
consists  of  two  steps. 

Stqp  1.  Deteimine  Icwer  bounds  (£,S)  and  \;pper  bounds 
(s’,s)  on  s*  and  S*. 

Step  2.  Find  the  oollectioi  of  all  (s,S)  policies  that 
minimizes  g(s,S  -  s)  -  over  the  cl^s  of  (s,S) 
policies  falling  within  the  bounds  defined  in 
Sti^  1.  Any  policy  in  the  collection  is  optimal. 


An  interval  bisection  or  a  Fibonacci  search  technique  is  then  enplcyed 
to  find  the  optimal  values  of  s*  and  Q*  =  S*  “  s*.  Note  that  this 
procedure  does  ncTt  guarantee  the  optimality  of  g(s,S  -  s)  since,  in 
general,  g(s,S  •>  s)  is  not  a  unimodal  function  in  s  and  Q. 


aa 


BtasaiAiiB&aai 


Ibr  other  work  done  in  the  area  of  one-oarnodity  inventory, 
the  survey  articles  of  Scarf  19]  and  Veinott  [18]  should  be  noted. 

Let  us  now  turn  to  the  ra-ocnnodity  (m  ^  1)  problem.  As  most 
practical  problems  in  inventory  involve  more  than  one  product,  there 
has  been  substantial  interest  recently  in  studying  the  theory  of 
rnulti-cxxtinodity  inventory  problens.  Ihe  problem  we  i^all  consider 
is  the  familiar  dynamic  inventory  system  with  periodic  review  [17] . 
Demands,  {D^},  for  the  items  over  a  sequence  of  successive  periods 
(i  *=  1,2,...)  are  assrtned  to  be  described  by  a  sequence  of  non¬ 
negative  randan  variables,  independently  and  identically  distributed. 

No  specific  restrictions  will  be  placed  vpon  the  statistical  properties 
of  the  demands  for  ocnmodities  within  a  given  period;  tiius,  in  general, 
the  danand  distribution  is  given  by  a  joint  density  function 
(t  Oatplete  backlogging  and  jjimediate  delivery  of  orders  are 

assunned,  and  as  previously  discussed  the  ordering  decision  in  each 
period  is  affected  by  a  single  set  tp  cost,  a  linear  ordering  cost  and 
holding  and  shortage  costs.  Follcwing  Sivazlian  [13] ,  we  define  the 
following: 

In  E?",  let  fi  =  {x|x  S)  vhere  x  =  (x^  ,X2 , . . .  ,3^)  is  the 
inventory  levels  of  all  items  prior  to  a  decision  and  S  =  (Si,S2,..wS^) 
is  a  non-negative  vector.  Let  the  lypersurface  that  subdivides  the 
set  fi  into  the  subsets  a  and  o®  be  defined  implicitly  by  an  equaticxi 
of  m  variables,  Z(S  -  x^)  =  0,  x^efl,  where  the  real-valued  function 
2(S  -  x)  has  the  following  properties: 


1.  Z(S  -  x)  is  defined  and  continuous  for  V  xefl 
2*  Z(S  -  x)  <  0  V  xea° 

3.  Z  (S  -  x)  >  0  V  xea 

Define  T  =  {xjxen,  Z(S  ~  x)  =  0}  and  U)(x)  =  x;  t,  x  eo®}; 

The  stationary  policy  considered  is  to  order  all  items  for 
saa®.  For  the  case  of  two  ocmmodities,  Pig.  (1.2)  illustrates 
geometrically  the  decision  regions.  It  is  inportant  to  note  that 
for  x^,  the  (a,S)  policy  is  not  feasible.  Hcwever,  once  a  stock 
level  readies  a  level  xeiJ,  the  policy  is  and  remains  feasible. 

With  this  ordering  policy,  the  sequence  of  inventory  levels 
Xi,X2,.../X  ,...,  vbere  X,  is  the  inventory  level  of  all  items  at 

•  JL 

the  beginning  of  the  i^h  period,  forms  a  Markov  process  since  the 
demands  are  independently  distributed  frcm  one  period  of  time  to  the 
next  and  since  the  (a,S)  policy  depends  only  on  the  last  state  of  .the 
process.  Sivazlian  [10]  and  [13]  studies  the  stationary  behavior  of 
the  system  for  the  special  case  viien  P  is  admissible,  i.e., 
rnA(x)=  4) (null  set)  vhere  A(x)  =  {tjx  <  t  <  S,  xeo®}.  If  W(y) 
denotes  the  distribution  function  of  the  stationary  distribution  of 
of  the  stock  levels  imrediately  after  a  decision,  then 


Pig,  (1.2).  (a,S)  ordering  policy  for  a  tvro-ccnnodity 

inventory  system. 
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Sivazlian  [10]  shews  that  the  solution  of  the  integral  equation  is 
given  by 

00 

V»(y)  =M  I  (S  -  y) 

n=l 

vhere  4''  ^  (x)  is  the  ordinary  n^h  fold  oanvolution  of  4)(x)  with 
itself  and 


_ 1 _ 

1  +  f  4>(t)  ^ 

*S“tew(S) 

Bie  optimal  decision  variables  Si,  S2,...,Sj^  and  the  (X)nfiguration 
of  the  adnissible  set  F  are  derived  by  minimizing  the  total  e3q«(7ted 
oost  per  period  given  as 


JU 

g(a,S)  =E(C)  =KM+  I  [c.E(D.)  +  h.E(y.)  +  (h.  +  p.)  E(B.)]  (1.3) 

••  i“l  1  i  1  i  XX  X 


vhere  C  is  the  total  enst  per  period;  and  for  items  i  =  l,2,...,m 
B  variable  cost/unit 

hji^  =  holding  oost/unit  at  the  end  of  a  period 
p^  B  shortage  oost/unit  at  the  end  of  a  period 
B  unfilled  demand  at  the  end  of  a  period 


B  inventory  level  at  the  end  of  a  period 
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£{•)  denotes  the  e^^ctation  operator.  > 

Fca:  the  tero-oonnodity  prbblan  when  the  denands  for  the  itene' 
are  independently  and  exponentially  distributed  with  '=  ■ 

0  tj^  <  “>f  Xi  >  0  (i  =  1,2) ,  'the  optimal  policy  is  investigated  for 
the  two  given  classes  of  adnissible  T 

;  *  I 


Case  1.  Here  F  is,  tlie  set  of  points  fir  vMch 

Xi  (Sj  -  xi).  +  X2(S2  -  X2)  =  a  ' 

v^ere  a  is  an  unkna^  positive  , 

I  *  , 

parameter  to  be  determined. 

!  ■ 

Case  2.  Here  F  is  the  set  of  points  fop  vMc*  ' 

max  IXjtSj.  -  Xj)  -  bj,  X2{S2  -  X2)  -  b2)  =  0  1 

viiere  bj  and  b2  are' unknown  j^sitive  ' 

I  I 

parameters  to  be  determined. 


IV-:  ’  Jtli  cases  -the  optimal  decision,  variables  are  ocrputed.  Moreqver 

J 

iJi^-azlian  shews  by  numerical  exanples  that  the  joint  (a,S)  policy 

I  I 

yields  a  ru  or  epprating  cost  than  the  case  vrfien  each  itan  is  ordered 

separately.  In  another  paper  Sivazlian' [12]  considers  the  multi- 

1 

oornodity  control  prcblera  operating  under  a  stationary  {o,S)  policy. 

•  .  '  *  ' 

Bie  joint  density  function  of  demand  for  all  items  in  a  given  period 

,  I 

?  I 

is  given  by 


m 

4»(t)  =  n  li*.  (t.)  0  <  t .  <  »,  (i  =  1,2, ...m) 

-  i=l  i  i .  -  -A  ' 


I 


I 


I  **A  *  * 

v^cre  4>.  (t.)  =  X.e  ^  ^  (i  =  l,2,...,m) .  F  is  adnissible  and  here  it; 
X  ^  ^  ‘ 

is  the  set  of  points  for  vMdi 


Xi (Sj  *  Xj)  +  X2(S2  “  X2)  +  ...  +  =  a  a  >  0 


In  I  this  case,  let 


=  [  r;>(t)  ^ 

Js-te<i){S) 


By  generalising  the  concept  of  Dirichlet's  multiple  integrals,  Sivazlian 

•  I 

shcK^  that  A  can  be  ej^ressed  as  a  single  integral.  Next/ similar 
e}?)ressions  for  E(Yj^)  and  E{B^)  are  obtained.  For  the  case  vhen  the 
variable  a'  takes  on  large  values,  ejqplicit  analytic  solutions  for  the 
Cptimal  values  of  S  and  a  are  obtained  by  minimizing  (1.3) ,  It  is 

•'  1 

eilso  shown  that  v^ien  m,  ^e  nuiber  of  itans,  becomes  very  large  and  a 
I  is  finite,  the  probability  of  placing  an  order  in  any  given  period  tends 
to  1,  and  tiiat  it  is  alaost  cptimal  to  use  the  "order  vp  to"  policy*, 
fbr  the  two-ocmmodit^'  problem,  Wheeler  [20]  follows  Iglehart's  [7] 
approach  for  the  one-ooranodity  (s,S)  model  in  showing  that 

I  ^  n 

vhere  is  the  total  e:q)ected  cost  of  operating  the  system  for  n 

periods,  given  that  x  is  the  initial  inventory  level  and  an  optimal 
policy  is  used  in  each  period,  and  g(o*,S*)  is  the  optimal  stationary 
total  expected  oo^t  per  period.  Follcwing  Greenberg's  [5]  notations 
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and  procjedure,  Wheeler  eilso  attenpts  to  find  the  tTcUisient  solutions 
for  the  probability  distribution  of  the  stock  levels. 

WacJcer^  fomrulates  and  studies  the  riiulti-ocntnodity  inventory 
problem  with  periodic  review  when  a  stationary  m-dunsnsional  (s  ,S) 
policy  is  followed.  Following  the  work  of  Roberts  [3]  he  gives 
analytical  fonnulas  for  determining  the  asynptotic  reordering  region 
viien  the  set  \:p  and  shortage  costs  have  large  values. 

For  the  case  viien  the  demand  D  in  each  period  is  discretely 
distributed  with  probability  mass  function  ,  viiere  x  is  the 

stock  level  at  the  beginning  of  the  period,  and  <|>0(j£»O)  <  1,  Johnson 
[8]  uses  the  policy  iteration  method  to  choose  the  optimal  (a,S) 
reorder  policy  vhich  minimizes  the  stationary  total  expected  cost  per 
period.  Further  he  provides  a  ccnputational  algorithm  for  finding 
bounds  and  diaracterizing  the  optimal  policy.  The  algorithm  is  a 
search  procedure  and  only  attenpts  to  find  the  configuration  of  the 
optimal  policy  pointwise. 

Finally,  the  work  of  Soland  [15]  is  briefly  mentioned.  In 
solving  seme  renewal  processes  in  two-dimcaisions,  Soland  solves  a . 
continuous  review  inventory  problem  involving  two  products  for  the 
case  vhen  the  demand  is  discretely  distributed.  In  minimizing  the 
total  expected  cost  per  cycle  (defined  as  the  time  between  ordering) 
he  concludes  that  the  optimal  policy  is  to  order  both  products  at  the 
end  of  a  cycle. 

Here  again,  as  in  the  one-ooimodiiy  problem,  a  renewal  type 
integral  equation  in  m-dimensions  is  solved  by  each  of  the  authors  in 
carrying  their  analysis. 


*See  additional  references 


A  nore  explicit  definition  of  the  problem  is  given  next. 

1.3  Definition  of  the  Problem  eind  £Ui  Outline  of  tlie- Chapters 

A  ntulti-oarniodity  inventory  system  vdth  periodic  review 
operating  under  a  stationary  (a,S)  policy  is  considered.  Ihe  ordering 
decision  in  each  period  is  affected  by  a  single  set  tp  cost  K,  a 
linear  variable  ordering  cost  c  =  (cj ,  C2» . . .  ,  and  the  ejected 

holding  and  shortage  cost  function,  L(x) ,  conditioncil  on  being  in 
stock  level  x  =  (xi,  X2,...,J^)  at  the  beginning  of  a  period;  L(x)  is 
assumed  to  be  tavice  differentiable.  Demand,  {D^},  for  the  items  over 
a  sequence  of  periods  (i  =  1,2,...) ,  is  assuned  to  be  independently  and 
identically  distributed  continuous  non-negative  randan  variables  with 
emtinuous  joint  density  function  <j>jj(t) .  Delivery  of  orders  is 
ijimediate  and  total  backlogging  of  unfilled  demands  is  assumed. 

Under  a  stationary  (o,S)  .  policy,  as  discussed  before,  either 
eiU  itens  will  be  ordered  to  bring  the  inventory  level  to  S  =  (Si,  Szf 
...,Sj^)  if  X,  the  inventory  level  at  the  beginning  of  a  period  prior 
to  ma3dng  a  decision,  is  in  a  (ordering  region)  or  nothing  is  ordered 
if  xeo®  (not  ordering  region) .  The  principal  concern  of  this  study  is 
to  find  the  cptimality  condition  for  (a,S)  policies.  Hds  is  done  by 
minimizing  the  ej^ression  for  the  stationary  total  expected  cost  per 
period  with  respect  to  the  decision  variables  that  characterize  the 
policy  being  used. 

3he  stutity  tegins  in  Qiapter  II  by  introducing  seme  mathematical 
concepts  .that  are  needed  in  subsequent  chapters.  First  the  concept 
and  properties  of  a  generalized  convolution  operator  defined  on  locally 
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integrablc  real-valucsd  functions  of  m  (m  >  1)  variables  are  introduced; 
next  the  solution  of  an  integral  equation  of  the  renewal  type  in  m 
.  dimensions  is  studied. 

Chester  III  the  analysis  for  deriving  an  analytical  e}q)ression 
for  the  stationary  total  eaq^ected  cost  per  period,  proceeds  in  three 
st^s:  (1)  the  basic  functiaial  equations  which  relate  the 

■  total  e:5)ected  cost  of  operating  the  system  for  n  periods  vhen  starting 
with  and  h^(x) ,  the  total  expected  cost  of  operating  the  system 

for  n  periods  ^rflen  starting  with  xea,  are  derived;  (2)  for  large  n  it 
is  shown  that  =  ng  +  u(x)  and  h^{x)  =  ng  +  v(x) ;  vhere  g  is  the 

stationary  total  e:5)ected  cost  per  period  and  u(x) ,  v(x)  are  functions 
of  the  initial  stock  levels;  and  (3)  an  analytical  expression  for 
g,  vMch  is  a  function  of  {a,S) ,  is  derived  by  solving  a  renewal 
type  integral  equation  that  results  from  the  use  of  the  asynptotic 
,  e>5jressions  for  fj^(x)  and  h^{x)  in  the  functioiial  equatioiis. 

In  Chester  IV,  the  optimality  conditions  for  {a,S)  policies 
*  are  studied.  Ihe  minimization  of  the  e35>ression  for  the  stationary 
total  ejqiected  cost  peir  period  with  re^cL  to  the  decision  variables 
proceeds  as  follcws;  First,  it  is  shewn  that  at  optimality 
r*  {xjxefl;  L(x)  -  C*  =  0},  vhere  C*  is  the  nunimum  value  of  g(a,S) 
excluding  the  variable  ordering  cost.  Ihen  the  necessary  and  suf¬ 
ficient  conditions  for  the  existence  of  proper  relative  minima  of 
9(0 ,S)  at  the  pair  (S’*',C*)  are  established. 

In  Chapter  V,  the  study  considers  the  two-ootmodity  inventory 
control  problem  viiere  the  demands  for  each  oonmodity,  in  a  given  period. 


are  independent  and  exponentially  distributed,  and  vbere  the  holding 
and  shortage  costs  for  each  oonno^ty  arc  linear.  The  equations  used 
to  solve  for  the  optimal  policy  parameters  are  analytically  determined. 
In  a  numerical  exanple  a  cniputer  program  is  developed  to  determine  . 
Si*,  S2*  and  C*. 

Finally,  in  Chanter  VI  recoiuiendations  for  future  research 
are  presented. 

As  oaipared  to  other  vrorks,  in  the  m-cmmodity  (m  >  1) 

-Ji 

inventory  control  problem  operating  under  a  stationary  (o,S)  policy^ 
this  research  is  more  general,  in  the  sense  that  no  assunptions  are 
pvade  about  the  ooxafiguration  of  the  ordering  region  or  the  admissibility 
of  r.  The  oily  restriction  an  the  set  T  =  {xjxen,  Z(S  -  x)  =0}  that 
subdivides  the  set  ft  into  the  ordering  rugicn,  0,  and  the  not  ordering 
region,  a°,  is  that  the  real  valued  function  Z  (S  -  x)  is  continuous 
and  defined  for  all  xeft.  Itoreover  the  ordering  r«^ions  must  be  siitply- 
cxjnnected.  The  ordering  regions  considered  by  other  researchers,  in 
the  m-oonroodity  problem  for  continuous  demand,  are  all  special  cases  of 
■the  general  case  considered  here.  In  E^,  Fig.  (1.3)  shows  gecmetrically 
the  configuration  of  seme  of  the  regions  that  are  considered  in  this 
research.  Fig.  (1.4)  shews  seme  ordering  regions  that  will  not  be 
considered.  The  m-dimensicmal  ordinary  convolution  operation  used  by 
other  researchers  in  stxxiying  the  ^stem  umder  consideration  is  a 
special  case  of  the  generalized  convolution  operation  introduced  in 
Chester  II.  For  the  case  when  the  set  F  is  admissible,  the  generalized 
convolution  of  any  function  reduces  to  the  ordinary  convolution  and 


. . 
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all  the  generalized  convolution  operation  properties  with  the  ccmnutative 
law  hold.  The  cinalysis  used  in  deriving  the  egression  for  the  sta¬ 
tionary  total  ej^jected  cost  per  period  is  shnilar  bo  Sivazlian's  [11] 
approach  for  the  ono-cxxTnodity  problem  \^en  operating  under  a  {s,S) 
policy.  The  asyupbotic  e;5)rossions  for  the  costs  of  operating  the 
system  are  deduced  hy  appealing  to  the  results  of  Yosida  and  Kakutan 
[22] .  nie  set  of  equations  that  will  be  used  to  determine  the  cjptimal 
parameters  of  the  adopted  policy  are  similar  in  form  to  the  set  of 
equations  used  hy  Sivazlian  [11]  in  determining  the  optimal  values  of 
s*  and  Q*  =  s*  -  s*  for  the  ont  jcranodity  problon. 


aiAPTER  II 


mTHEMATICTVL  PRELIMINARIES 

2»1  Introduc±ion 

In  this  chapter  a  number  of  mathematical  conc^ts  that  will 
be  needed  in  subsequent  ch^ters  are  discussed.  First#  the  concept 
and  properties  of  a  generalized  convolution  operator  defined  on 
locally  integrable  real-valued  functions  in  m(m  ^  1)  variables  are 
studied  and  then  the  solution  of  a  renewal  type  integral  equation  is 
given. 

2.2  Ihe  Generalized  Convolution 

Let  S  be  a  positive  vector  in  a  ^aoe  of  m{m  ^  1)  dimensions; 

i.e.#  S  >  £  <“=>  >  0/  i  =  l#2#...#m.  In  define  the  set  ft  = 

{x  I  X  S}  #  and  let  the  hypersurface  that  subdivides  the  set  ft  into 
the  subsets  a  and  o°  be  defined  inplicitly  by  an  equation  of  m 
variables;  Z(S-x°)  =  0#  x°eft#  where  the  real-valued  function  Z  (S-x) 
have  the  foUcfwing  properties; 

1.  Z  (S-x)  is  defined  and  continuous  for  V^ft; 

2.  Z (S-x)  <  0  VxGO® 

3.  Z  (S-x)  >  0  Vxea 

Next#  define  the  sets:  (a)  w(S)  «  {x  |  JKft;  Z  (S-x)  <  0}; 

(b)  w(x)  s*  {t  I  t#  xeu(s);  t£x}; 

(c)  r  =  {x  I  3«ft;  Z(S2X)  =  0); 

(d)  w'(t,x)  =  {v  I  £<v<x;  t,vew(x)}; 
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(e)  r^Cx)  =  {t  I  t  <^x;  }few(S)); 

(f)  r'(x)  =r«(x)  -w{x) 

Let  us  make  the  transfonnation  u  =  x-ty  tew  (x) ,  then  the 
set  R(x)  is  the  image  of  the  set  w(x) ,  is  the  image  of  F,  and 
Rj  (x~t)  is  the  image  of  w'  (tyx) .  In  the  particular  case  m  =  2, 

Pig,  (2.1)  illustrates  geometrically  the  effects  of  the  transfor¬ 
mation  on  the  above  defined  sets. 

For  ucR(x) ,  define  the  set  R^  (u)  =  (t  |  0  ^  t  <_  u;  teR(x)}. 
Consider  the  set  of  functions  D,  f(t)eD  is  a  real  valued  function  of 


m  variables  (t,  yt,  y . . .  yt  )  for  vMch  we  have 
*  m 


1,  f  (t)  a  0  for  t  <  ^ 

2.  f  (t)  is  locally  integrable  on  Rj  (u) 


(2.1) 


Cn  this  set  of  functions  y  addition  and  g-convolution  operations 
are  defined  as  follows: 


.Defx 


1.  Addition  (f4g)  (iy  *=  f  (j^  +  g(t); 

2.  g-Convolution  (f*g)  (u)^^  [f  (u-v)  g(y)  dy 

ieR7@ 


then  the  following  relations  are  satisfied 

1.  f(t)y  g(t)eD  =>  f(t)  +  g(t)eD 

2.  f  (t)  +  g(t)  *=  g(iy  +  f  (t)  (Cornutative  Law) 

3.  {f(iy  +  g(iy}  +  h(t)  “  f(t)  +  {g(t)  +  h(t)}  (Associative  Law) 

where  the  equality  sign  is  used  in  the  ordinary  sense.  And  the  set 

of  functions  D  is  an  abelian  group  vd.th  respect  to  addition  in  which 
the  zero  element  is  the  function  f  (^)  s  0  and  the  inverse  of  the 
function  f (t)  is  -  f (t) . 


3 

I 


•.k-<  r.  .. 


2 


Properti 


of  q-Convolution 


1.  ■  Closure  Lav; 

If  f(t),  g{t)d5#  then  h(u),  given  by 


h(u)  »  [  f  (uH:)  g(t)  ^ 


belongs  to  D. 

Proof  t  h(u)  is  locally  integrable  on  R|(u)  since  both 
functions  f  and  g  are  locally  integrable  on  Rj  (u) . 

Moreover  h(u)  is  idoitically  zero  for  u  <  0^. 

2.  Distributive  Law 

If  £,  g,  h,  e  D  then 

If  *  (g+h)l  (u)  =  (f  *  g)  (u)  +  (f  *  h)  (u) 

Rjoof:  definition 

If  *  (g  +  h)l  (u)  =  f  f(u-y)  [g(y)  +  h(^)J  ^ 

•  1  f(u-y)  9(Z>  ^ 

^ycR,  (u) 

+f  f  (u-y)  h(y)  ^ 

'  XeR,(u) 

«*  (f  *  g)  (u)  +  (f  *  h)  (u) 

3.  Associative  law 

If  f ,  g,  c  D  then  [f  *  (g  *  h)l  (u)  «  I  (f  *  g)  *  hi  (u) 

Proof;  By  definition 

If  *  (g  *  h)l  (u)  «  [  f  (u-v)  I  g(v>-y)  h(y)  dy  dv 

(y) 

•  f  ff  g(v-y)  h(y) 


(P-2 


vjhere  { (v,y)  |  VeRj  (u) ,  (v) )  which  Is  equivalent 

to  the  set  defined  as 

^  I  Z  <  3L  <  H.?  V,  XpRj  (u)  > 

Also  by  definition  we  can  write 

I(f  *  g)  *  hHu)  =  I  I  I  f  (u-yt)q(t)  h(j;)  ^ 

(H)  teR^  {u~y) 

Let  t  ®  V  “  then  (2.3)  beoones 

((f  *  g)  *  h]  (u)  =  I  (  I  f(u~v)  g(v-y)  h(jO  ^ 

Rj  (u)  v^Rj  (u-y) 

Since  g(v-y)  =  0  for  0,  v-y:R^  (u-v)  iirolies  that  jjrR*  (^,u) , 
where  «  (v  |  X.1.  Zl.  E»  *  ^ence  we  have 

((f  *  g)  *  hl(u)  =1  I  I  f(u--v)  gCv-y)  h(^)  ^ 

«|  I  f (u-^)  g(v~y)  h{^)  dv^ 

W, 

>Aiere  W ,  =  ftv»Z)  I  Zl.  Z  -  Z'  Z»2§^  •  In  E  *, 

Fig.  (2.2)  illustrates  geonvetrically  the  sets  appearing 
in  (2.2),  (2.3),  and  (2.4). 

Ihe  sets  W ,  and  W  j  as  defined  are  equivalent,  and  that  in 
turn  Inplies  W  ^  and  W ,  are  equivalent.  Since  the  integrals  given 
hy  (2.2)  and  (2.4)  exist,  then  by  Pubini's  theorem  the  iterated 
integrals  are  equal.  Iherefore,  for  eadi  ueR(§) 
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4.  C3otmiutative  Lav;: 

■{-»■■'■■■■» . . 

We  observe  that  in  general,  if  f ,  geD,  then 
(f*g)(u)  (g*f)  (u) 

definition 

(f*g)  (u)  =  I’f  (u  -  t)  g(t)  ^ 

,  teRi  (u) 

Let  u  t  =  then  : 

(f*g)  (u)  ^  I  f  (V)  gM  (g*f)(u) 

veRi  (u' 

viiere  Ri(u)  =  {v|u-veRi  (u) ). 

,  Note:  ■  the  g-oonWlntion  of  f eD  with  itself  cxamniLiS. 
Moreover,  if  we  define  for  n  =  1,2,,.. 

;  f(l)  (u)  =  f  (u)  , 

=  j  *to-l)<H±>  ‘4’  «  (">2) 

teRi  (u) 

I 

then  from  the  g~oonvolution  associativity  law  we  can  write 

,  j 

teRj  (u)  tcRi  (u)  (n  ^  2) 

If  we  let  u-t  =  y;  then  Rj  (u)  is  mapped  into  Rj  (u)  and  we 
■  can  v^ite  '  *”  "* 

“  I  ^(n-1)^— ^  f  ^(n-l)^~  — 


•  i.eRi(u) 


tcRi  (u) 


:  |  fM  ^=1  f(t)  6t 


tcRi  (u) 


tcRi  (u) 


(n  1  2) 


If  Pig.'  (2.3)  illustrates  geanstrically  the  sets 
R^  (u)  and  Rj  (u) . 


WjrWWrVe''* 


J'l” 
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Ifenceforth  we  conclude  that  the  set  of  fiincitions  D 
defined  by  (2.1)  forms  a  non-comtutative  ring  with 
respect  to  addition  and  convolution  ot^erations. 

Remark  2.1; 

It  is  important  to  note  that,  in  the  case  vdien  for  all 
ueR(S) ,  Ri  (u)  =  {t  I  ^  £  t  5. the  ordinary  convolution  in  in  (m  ^  1) 
dimensions  is  defined  and  all  the  previous  properties  with  the 
Ccnnrutative  Law  hold.  Thus,  in  a  sense,  the  ordinary  convolution  is 
a  special  case  of  the  g-convolution.  It  may  also  be  noted  that  for 
the  case  v^n  Tj  is  admissible  the  g-con\»olution  of  any  function 
reduces  to  the  ordinary  convolution. 

Remark  2.2; 

Let  (j>(x)  be  a  locally  integrable  density  function  of  the 
irandora  variable  D.  For  every  xew(S) ,  define 

♦  (DW  =  I  (»)(iy  ^ 
x~teu)(x) 

X 

=  p(t)  ^ 

and  for  n  ^  2  . 

”  *(n)^2l)  =  J  ^(r.-i)(x-t)  6t 

x-te(i)(x) 

X 

«!"^(x^)  =  I  *(n-l)(^)  ^(t)  ^ 


Pmn  the  above  definitions  and  the  definition  of  the  set  w(x) ,  we 
have  ly  induction 


^(n)  W  1  M 


in  >  1) 


Ih’js, 


I  *(n)  1  ?  (x) 

tt=l  n=l 


By  definitiai 


t(x)  =  Ht)  . . *1—  *m 

r*i 

=  1^  <<>i(t^)  %  =  ^^(Xi)  (i=l,2,...,m) 

f 

where  «>i(xi)  is  the  marginal  distribuUon  of  the  randcm  variable  D^ 
with  density  function  4)i(xi) .  If  we  denote  (x^)  the  n-fold 

ordinary  convolution  of  ,  then  for  n  =  1,2,... 

(x)  <  $.  (X.) 

^  — . . .  »  ♦  • 

•ft* 

4*^2)  (x)  -  f-  $<i)(x-t)  <j)(t)  dt 

Jo  - 


i  ^ 


4  ^ 

-J 


Hie  integrals  given  by  (2.8)  and  (2.9)  exist  .  Hierefore,  by  . 


i  I  : 

Ftibini's  theorem, the  iterated  integrals  are  equal.  Hence  for  xea)(S) 


v>e  have 


a(x)  =  b(x) 


From  (2.6)  since  ^(v)  =,  0  foif  v^R,  (x-y) ea)‘(x-y) ,  we  have' 


b(x)  -  \  .  if  'f(x-y-v)'g(v),dv  1  h(y)  ^*C(x) 


And  this  cqnpletes  the  proof  of  the  leirma. 


Note?  Ifit  x-y-v  =  t,  then  fron  (2 .,7) 


<1  ' 

C'(x)  =  f  ■  .  (  f  .  f (t)  '4(x~y“t)  ^  ’)  h(y)  ^ 

•'yea)(x)  ^teR,  (x~y) , 


Remark  2.3;  let  ((.(t)  be'  a  locally  integrable  density  function  of  the 

■  '  '  I  ‘ 

randan  variable  Dr  and  let  4,  .  (t),  n  >  1,  be  the  n-fold  g-oonvoluticai 
of  <>(^)  with  itself.  From  the  g-oonvolption  properties  we  have 


,  f  (J>(x-y-v),  4,  .  (V)  ^ 

W,(if2)  - 'veRiCJoS 


'♦'rnito'  (”>1> 


Hence  from  l£imti  2.1 


'  i 

f  .  (^(x:;t),  I  f  h(y)  dt 

^teu)(x)  (5^  i) 


'yGw(t) 


r  ,  ♦(n)<2ac>  h(z)  a: 


J  .1 

I 
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Similarly  from  the  g-cx)nvolutic«  properties  and  Lamia  2.1  we  have 

e 

f  f  h(^)  ^  dt 

JtewU)  ^yea3(t)  — 

V  < 

=  1  ,  *(n)^2SrJi)  h(ji)  ^ 

•'yea)(x)  '  ' 

In  E*,  Fig.  (2.4)  illustrates  geometrically  the  sets  a)(x) ,  to(x-y) , 
and  Ri  (x-y) . 

Theorem  2.1;  Let  (j)(x)  be  a  locally  integrable  density  function  of 
the  random  variable  D.  Then,  for  every  xetu(S)  and  |1a(x)  ]  <  «>  there 
exists  a  lanique  solution  to 

U(x)  =A(x)  +  f  <j)(x:;y  U(t)  dt  (2.10) 

*“  'teaj(x) 


given  by 


U(x)  =  A(x)  +  ij;(3^)  A(y  ^  (2.11) 

'tea)(x)  ^ 


vhere  ij;(x-t)  =  I  (x-t) ,  (x-t)  is  the  n-fold  g-convolution  of 

1V"1 

<|>(x-t)  with  itself,  and  |^(x)i|  =  Sup  |a(x)  |. 

»EW(S) 

Proof;  To  prove  solution  existence  and  uniqueness,  successive 
approximation  method  is  used.  Consider  the  sequence  of  functions 
defined  inductively  as  follows: 


Uo(x)  =  A(x) 

Ui(x)  =A(x)  +  I  Uo(:W  ()>(x;:l)  ^ 
^tEW(x)  " 


and  U  (x)  =  A(x)  +  f  A (x-t)  ^  1  0) 

'te(;j(x) 

^^lere  the  subscripts  refer  to  the  inductive  process.  Let  Ajdenotes 

the  collection  of  all  such  functicais.  Aj  forms  a  Banach  space  in 

which  the  norm  of  an  element  is  defined  by 

■  ||U(x)||  =Sup  lU(x)| 

xeu(S) 

To  proceed  in  the  proof,  first  v,’e  wish  to  shov*;  that  the  sequence  of 
functions  converges  to  a  function  U(x)  for  all  xew(S).  Next, 

we  shall  show  that  this  function  U(x)  satisfies  (2.10),  and  finally 
that  U(x)  is  a  unique  solution  to  (2.10). 

Fran  the  recurretice  relations  we  can  write 

U  (x)  -  U  (X)  =  [  (Kx-t)  [  U  (t)  -  U  ,  (t)  )  ^  (n  >  1) 
n+l-  n-  -  n-  n-1  -  -- 

Ta)dLng  absolute  value  of  both  sides  we  obtain 

|Un+lW  "  ^n^  I  ^  f  |Un<t>  ”  Vl^^)  1  ^  (n  >  1) 

'tew(x) 

Let  m  =  |A(x)  | 

Upon  taking  norms  of  both  sides,  of  the  above  relations,  we  get  for 

n  ^  0,  1,  2,  .  •  • 

|U.  (x)  -  U  (x)  I  <  mf  (^(t)  dt  =  m*.  .  (x) 

“  ^x-tefa)(x)  “ 
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lb  show  that  (2.11)  satisfies  (2.10)  we  proceed  as  follows: 

Using  (2.11)  we  can  write 

UW  -  [  6(3c^)  U(t)  ^=A(x)  +  [  il>(x"t)  A(t)  ^ 

Jt^(x)  *“  'tea)(x)  *” 

-  f  (Kx-t)  (A(t)  +  f  i|»(t:l,)  A(l^)  d^l  dt 

Jtp«(x)  -  Jt,ca(a  — *  - 

Applying  the  g-convolution  properties*  Lemna  2.1*  Remark  3.3*  and  the 
definition  of  ^(t)*  we  get 

0(x)  -  f  ^(x-t)  U(t)  dt  =  A(x)  +  y  [  .  (x-t)  A(t)  dt 

"  ■»teu(x)  ”  1^1  Jtew(x) 

”  [  6(x~t)  A(t)  ^  +  f  ^(JS^t)  • 

'^(d(x)  “  -'^(x) 


iff  A(^)  d^)  dt 

nil  't,6(o(t) 

«*  A(x)  +  y  f  .  (x-t)  A(t)  dt  -  f  (^(x-t)  A(t)  dt 

“  n=l  ^tew(3c)  '^(x)  ~ 

+  I  f  A(t)  ^ 

nSl  W(x) 


*  A(x) 

Thus  (2.11)  s\tisfles  (2.10). 

Conversely*  we  can  shew  that  if  U(x)  satisfies  (2.10)*  then 
0^)  satisfies  (2.11).  Using  (2.10)  we  can  virite 


A(Jc)  ♦  J  i(>(x~t)  A(iyi  ^a'U(x)  -  I  ♦(x-t)  U(t)  |dt  i 
'^(x)  ^^(x) 


+  f  {  U(t)  -  f  ♦(t-t;)  U(t,)  dt,]  dt 


Using,  again,  Lenina  2.1,  Ranark=2.3,  and  the  drfinition  of  we  rjet 
A(x) 


(x)  +  f  ii/Cjrt)  A(t)  dt  «  U(x)  -  f  (Kx-t)  rj; 
'tew(x)  -  ^tew(x) 


€0 

I' 


xiFl  ^tew(x) 


♦  (n)^25=t)  U(t)  ^  ) 


•  r 
- 1 1 


npl  *tcu(x)f 


U(t)  dt 


»U(x)  1 

.  ■  i 

ISius  the  solution  for  (2.10)  is  given  (2.11). 

New  if  the  set  F,  as  defined  in  Section  2.1,  is  admissible, 

*  I 

*^(n)  ~  v^iere  (iy  is  the  n-fold,  n  ^  1,  ordinary 

convolution  of  (Ji(tJ  with  itself,  Bias  the  solution  for  (2.10)  will  be 
given  by 


I 


•  *  ' 


■'  Vvr  ■»V^ 


r:::  i7-Ti-'^^’^?rj.-t7«_?^  r'r-v-’<^ 


t  ! 
} 
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2,4  Other  Mathenatical  Concepts 


As  ve  slvdl  see  later  the  following  theorem  fzon  calculus, 
[23] ,  will  be  needed  in  Chapters  IV  and  V. 

Hheorqn  2.2;  Mean  Value  Iheorem  for  multiple  integreds. 

If  f  (t)  is  a  Qontinuous  function  on  a  corpact  set  T  for 


which  A{T)  *  I  dt  0,  then  there  exists  a  point  v  in  the  i 

!  »*rf”  ” 


K 


I 


I  of  T  such  that 

I  .  f  f  (t)  dt ‘=  f  (v)  A(T) 

Jm  T 


! 


Remark  2.4: 


From  the  theory  of  Linear  Operators  in  Banach  Space  the 
fii'Uoiong  results  will  be  needed  in  Chapter  III. 

'  i 

Define  the  operation 

I  (Ta)  (x)  =  f.  a(x-t)  (J)(t)  dt  +f  a(x)  ^ 

""  'x-tew(x)  'x-te  r*  (x) 


=  :f  a(x,t)  (|>{t)  ^ 

Jx-ter''(x)  “ 


5  [  <>(x--v)  a(x,v)  dv 

^ver"  (x)  “ 


xew(S) 


• 


With  =  1,  l”  is  the  iteraj»  of  T,  for  n  «  1,  2| 

-T  is  defined  for  all  functions  a(v)  vhidi  are  bcunded  for  all  ve3t"(x) , 


xe(i)(S).  The  collection  of  all  such  functions  forms  a  Ban^  (M)  Space 


in  which  the  norm  of  an  element  is  defined  to  be 


1  ^ 


i  J 


i  i 

3  a 

i  % 
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aip  |a(v)l 
yer"(S)  “ 

Fcoc  such  integral  ^leraticm  T  with  bounded  density  kemelr  Yosida  and 
Kakutan  [22]  shewed  that 

7^  «  5  T .  +  P” 

•*  B  +  n  •  1/  2f  •  •  •  • 

where  B  .is  the  ergodic  part  of  the  operation  and  P^^  is  the  non-ergodic 
part  of  7^#  k  is  the  number  of  proper  values  of  T,  (i  **  1,  2,  .  .  .»  k) 
is  a  conpletely  continuous  linear  c^jeration  and  P^  (n  ^  1)  is  a  conpletely 
continuous  linear  operation  (vMch  might  vanish)  with 


where  r  and  e  are  positive  constants  independent  of  n. 


CHAPHJR  III 

MATHEMATICM.  PORMUIATICN 

3,1  Introduc±ion 

An  m-connodity  (m  ^  1)  inventory  control  system  operating 
under  a  (o»S)  policy  is  considered.  At  the  beginning  of  each  period 
a  decision  to  order  or  not  to  order  is  node  depending  on  the  stock 
level  X  =  (Xj ,  Xj ,  .  .  • »  •  B'e  ordering  decision  in  each  period 

is  affected,  as  discussed  previously  in  Chapter  I,  by  a  single  fixed 
set-up  vT'st  K,  a  linear  variable  orderly  cost  c  »  (Cj,  Cj,  .  .  c^), 

and  the  expected  holding  and  shortage  cost  L(x)  conditional  on  being 
in  stock  level  x  at  the  beginning  of  a  period.  Datand,  {D^},  for  the 
items  ever  a  sequence  of  periods  (i  =*  1,  2,  ...  .),  is  assumed  to  be 
independently  and  identically  distributed  continuous  random  variables 
vith  joint  density  function  t  t  =  (tj,  tj,  .  .  •» 

Ibinediate  delivery  of  orders  and  ccnplete  backlogging  of  unfilled  demands 
are  assixned  for  all  conmodities. 

As  discussed  in  Chapter  I,  under  the  adopted  stationary  (a,S) 
policy,  the  sequence  of  stock  levels  at  the  beginning  of  each  period 
farms  a  discrete-time  Markov  process. 

Ihis  chapter  is  devoted  to  finding  an  esqsression  for  the 
stationary  total  expected  cost  p^  period  based  on  the  assumed  ordering 
policy. 
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-  ■  - - - - - 


periods  v^en  an  order  is  placed  initially,  i.e., 
xco. 


If  no  order  is  placed  initially,  i.e.,  xeo°,  then 


f^(x)  satisfies  the  functional  equation 


f »  W 


(n  =  l) 


f_(x)  =  T,(x)  +  [  f  ,  (x-t)  (^(t)  dt 

”  “  -  •’x-tewCx)  ”  ^ - - 


+  f  h  -i  (x-t)  (^(t)  dt  (n  >  2) 

Vter»(x)^^ - - 


However,  if  an  order  is  placed  initially,  i.e.,  xea,  then  i^(x)  satis¬ 
fies  the  functioned  equation 
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■■\  -r",  ""c*  '=t^4v>aV*; ►>*^.^''^7V'''‘*T^s-' 


To  proceed  in  the  proof,  first,  we  shall  siirplify  (3.7)  by 
substituting  for  a(x)  as  given  by  (3.5).  Next,  we  shall  find  an  ana¬ 
lytical  expression  for  ^^^(S)  and  use  it  in  finding  the  expression  for 
fjj(x) ,  SCO®,  for  large  n.  Finally,  we  shall  determine  the  expression 
for  hjj(x) ,  JKO,  for  large  values  of  n. 

Let  =«  (pj  ,p2  ,  .  .  .  ,p  ^) ,  vAiere  p^  is  the  expected  value  of 
the  random  variable  i  =  1,  2,  .  .  .,  mi  Ercm  (3.5)  and  the  defini¬ 
tion  of  the  sets  (d(x)  and  r*  (x) ,  we  have 

a(x)  =  f  (K  +  c*^(S-(x-t))]  (^(t)  ^ 

^x-ter*  (x) 

=  rCK  +  c'^(S-(x-t))3  <>(t)  dt 
^0  "■  ■" 

-  J  [K  +  c^(S-t))  (}>(x-t)  ^ 

^t€w(x) 

=»  K  +  c*^(S-(x-p))  -  I  (K  +  c*^ (S“t)  1  (x-t)  dt 

^tea)(x) 

Upon  using  the  g-axwolution  properties  and  Lama  2.1  we  ctotain  for 

i  ^  1,  2,  ... 

f  .  (x-t)  a(t)  f  .  (x-t)  {K  +  c^(S-(t-p)) 

.  'tew(x)  ^tew(x)  *“ 

-  [  (K  +  c'^(S-ti))  A(t-ti)  dti)  ^ 

'tie(i)(t) 


r 


ii  ! 

s  t 


I  i 

I  i 


I  { 

U 


I  i 


^  ; 


K 
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1  (t>  (x-t)  [K  +  C*’(S-  (^) )  ]  dt 

^te(o(x)  -  — 


‘  f  [K  +  c'^{S-t)l  dt 

^teu(x)  '•  ■“ 


and, 


f  6  (x-t)  a(t)  dt  =  f  ^  (x-^:)  [K  +  c^(S-(tnj))l  dt 

'tew(x)  (i)  (i) 


tew(x)  (i) 


'tea){x) 

T^xjn  using  the  above  relations,  we  obtain 
n-2  f 

a(x)  +  T  (t>...{x"t)  a{t) 

*■  i=l  ^tew(x)  “ 


“  f  IK  +  c'’’(S-t)l  dt  (i  >  1) 


dt  =  K  +  c'^(S-x)  +  c^v 


-  f  {K  +  c^(S~t) }  (i>(x-t)  ^ 

teuCx) 

n”2  .  m 

+  y  [K  +  c  (S-t+u)l  dt 

ih  JtGu(x) 

-  I  f  [K  +  c^(S-t))  dt 

^teu(x)  “J 

On  sinpliJfying  and  regrouping  terms  we  get 
.  n-2  , 

a(x)  +  y  (j),,.(x-t)  a(t)  dt 

”  iSl  'teu(x)  - 


»  K  +  c 


p  n— 2  , 

ttti+  I  I 

i»l 


(x-t)  dt]  +  c  (S-x) 

tcu(x) 


“  f  ♦(n-l)  ^25::^^  [K  +  c^(SHy]  ^ 

^teto(x) 

VSCM,  we  shall  proceed  to  find  an  analytical  expression  for 
^(S),  n  >  2. 

Erom  (3.4) ,  far  jpS,  we  have 
f  j  (S)  =  L(S) 

f„(S)  «=  L(S)  +  f  f  ,{S-t)  (J)(t)  dt  : 

"  -  -  WaoCS)  ”-l  — - 

+  f  tK  +  c’y  (Kt)  at 

^s-ter^{S) 


+  f  ,  (S)  (^(jy  dt  (n  >  2) 

'S^tcX*  (S) 


[  [K  +  c^l  (j)(t)  dt 

^S-ter* (S) 

T(f  .(t)){s)  =  f  fn-i^§::l>  ^ 

n-1  -  -  i^e„(S)  "  ^ 

+  f  -(s)  f  (Kt)  dt  (n  >  2) 

^S~ter* (S)  “ 

Note  T  as  defined  is  a  bounded  linear  c^seration.  Substituting  (3.10) 
and  (3.11)  in  (3.9)  v?e  get 

fj(S)  «  L(S) 

fpjS)  «  L(S)  +  T(f„_j^(t))  (S)  +  c(S)  (n  >  2) 


Let  us  denote  by 
c(S)  = 


. . .  -■■‘•i**jgieg 


(3.10) 


(3.11) 

I 


] 

•i 


(3.12) 


i 

j 

i 


i 

/! 


Fran  (3.12),  by  inducticn  fco:  n  »  2,  3,  .  .  we  have 
fjCS)  =  L{S)  +  c(3)  +  T(L(t))  (S) 
f,(S)  =  L(S)  +  c(S)  +  T(L(t)  +  c{t)  +  T(L(t)))  (S) 


and, 

n-2  . 

f-(S)  =  L(S)  +  c(S)  +  I  TTCLCt)  +  c{t))(S) 
n-  -  “  i=l  ”  ~  “  , 

+  a^"^(L(iy)(s)  (n>2) 

Fran  Yosida  and  Kakutan  [22]  we  have 

•r  =  I  T .  + 

j=i  3 

*=  B  +  p^  (i  «  1,  2,  .  .  .) 

where  B  is  independent  of  i.  Ihe  linear  operations  Tj(j  «  2,  .  . 

and  are  e5q)licitly  defined  in  Remark  2.4.  On  substituting  (3.14) 
in  (3.13),  we  get 

f.(S)  a  L(S)  +  c(S)  (.(BfP*)  +  (B+P*)  +  .  .  . 

n  —  —  — 

+  (BfP^^)l  L(t)+c(jy)  (S) 

+  lBfP”"^l  (L(t))(S)  (n^2) 


Upon  regrouping  terms,  we  obtain 
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f^(S)  =  L(S)  +  c(S)  +  (n-1)  B(L{t))(S) 


+  {n-2)  B(c(t))(S)  +  [  I  pi  1  (L(t))(S) 

i=l 


+  t  "f  1  c(t))  (S) 
i=l 


(n  >  2) 


(3.15) 


Rran  (3.6)  and  the  definition  of  the  sets  ii)(S)  and  (S) ,  vje  can  write 

b  ,(S,x)  =  f  ,(S)  f  (Kt)  ^ 

n  1  —  n  1  *x“tEr*  (x) 

«  f„-l(S)  (f“<Kt)  dt  -  f  ♦(x-t)  dt  1 
"  “  h  ■*teo)(x) 

=  f„-l(S)  (1  -  [  ♦(x-t)  dt  1  (n  >  2) 

"  ■>tew(x) 


and  for  n  >  2 


b  .(S,x)  +  i  [  'J>„^(2^)  b  (S,t)  dt 

^  i=i  Jteaj(x)  'i  ""i'i 

t  r  "'2 

«  Vi^s)  i\  -  J  <Kx-t)  +  I  f  ,  .  (s)  • 

^  "  ''tGW(x)  |_  i=l 

I  -  f  ^ 

^tew(x)  '  '  '^ea)(t)  J 

Cb  applying  the  g-convolutioi  properties  and  Lenina  2.1r  we  can  write 

n“'2 

Vi<£'2.>  +  M  ^  1  S. 

i»l  ^tew(x) 


i 


+  L(S)  +  c(S^)  +  B(L(t))  (S).  +  P*  (L(t))  fS)l 


+  L(S)  [B^"^  - 


iL(S^)  +  C(S)]  [P*  -  P‘  +  F*  -  P*  +  .  .  .  + 


+  -  B^^]  +  L(S)  [B*^2  _  +  B(L(t))  (S)  ■ 


.[  (n-2)  CP®  -  f»]  +  tn-3)  [P»  -  i'*J  +  .  .  . 


+  2tP”“^  -  P”"^]  +[F^“^  -  B«“2]  J  +  B(c(t))  (S) 


[(n-3)  (F®  -  P»]+  (n-4)  tP^  -?*]+...+ 

n**2 

-  p”"^]  ]  +  [  Z  p3j  (L(t))  (S)  IP®  -  F*] 
j=l 

n-3  .  ,  2  j 

+  I  2  P^l  (L(t))(S)  [F*  -  F®J  +  .  .  .  +  I  Z  P^l 
j=l  -  -  ,  j=i 

(L(t))  (S)  IB^'^  -  (S)  IF"-3  - 

n-3  j  n-4  j  . 

+  I  Z  P^l  (c(t))(S)  [P®  -F*]  +  I  Z  PJ]  (c(t))(S) 

j=l  "  j=l 

■ 

[pl  -  p2j  +  .  ,  .  +  pi  (c(t))  (S)  Ib”'"^  -  p”“^] ' 


I 


On  siniplifying,  we  obtain 


n-2  f 

b  ,  (S,x)  +  T  (}»...  (x-t)  b  ,  .  {S,t)  dt 

n-1  -  -  Jtea)(x)  - 

I 

*  L(^)  +  c(S)  r  c(S)  F*'”^  -  LCS)!^"^  +  (n-2)  B(L(t))  (S) 
n-2  . 

-B(L(t))(S)  I  +  (n-3)  B{c(t))(S) 
i=l  “  “ 

I 

!  n-3  .  n-2  . 

-B(c{t)){S)  X  P^+t  y  P^l  (L(t))(S) 

“  “  i=l  .  j=l  ■"  “ 

+  (  -P^^(L(t))(S)P*  -  r^“^(L(t))p2-.  .  . 


-  P*(L(t))(S)p”“^]  +  t  -P*'~3(c(t))(S)F'  -  P^"^(c(t)(S)P^ 


,  n-3  . 

-  .  .  ,  -  P»(c(t))(S)P"”'’]  +  1  X  P^]  (c(t))(S) 

-  -  j»l  ”  ~ 

I 

*  •>  T 

L(S)  +  c(S)  -  c(S)f’^~^  -  L(S)P^  +  (n-2)  B(L(t))  (S) 


-B(L(t))(S)  I  P^  +  (n-3)  B(c(t))(S) 

;  ”  “  i«=l  - 

n-3  ,  n-2  ,  ’  _ 

-B(c(t))(S)  I  P^+C  X  P^l  (L(t))(S) 
“*"1=1  i=l  "■ 

+  I  I  P^)  (c(t))(S)  -  I  p""^”^(L(t))(S)F^ 
i=l  -  i-i 

n-3  _  .  j 

-  y  (c(t))(S)F^ 

1=1  *"  “ 


(3.19) 
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'A'  ,  r'^^'T*.  ^  'r''»-'^’j  »'"’s^^'’^  77"'’'.'yi'i7’«i~  j ;  . . 


V2<i'2^> 


tir2  .  n-3  . 

[  I  P^]  (L(t))(S)  +  t  X  P^]  (c(t))(S) 
1=1  *"  “  i=l  -  ~ 


^  r^^“^(L(t))  (S)pi 

i=l  ~  “ 


“  I  p^^"^(c(t))  (S)F^  {n>2\ 

isl  ~ 


(3.20) 


Then  (3.19)  becxxnes  (after  substituting  for  ^1.1) 

n-2  / 

*=U-,<S,x)  +  I  <^/ix(x-t)  b  .  .(S,t)  dt 

nr-i - (i)  -  n-l-i - 

«  L(S)  +  c(S)  -  c(S)  J  (x-t)  dt 

"  ”  ~  ‘'tew(x)  '*'"2) 

-  L(S)  f  (x-t)  dt  +  (n-2)  B(L(t))  (S) 

”  ''tew(x) 

n**2 

+  (n-3)  B(c(t))(S)  -  B(L(t))(S)  I  f  (x-t)  dt 

-  ~  “  i«l  'teu(x) 

n-3  f 

••  B(c(t))  (S)  I  I  (x-t)  dt  +  D„  -(S,x) 

”  ■"  i=l  'tEw(x)  (i> - n  2  ~  - 


(n  >  2) 


(3.21) 


Using  (3.8)  and  (3.21)  in  (3.7)  we  get 


f„(x)  «  lL(x)  +  I  f 
"  "■  i»l 


i*»l  'tea)(x) 


♦  /.v(x-t)  L(t)  dtl 
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T  f 

+  [K  +  c  V  [1  +  I  <l>/i^{x-t)  dt  1 

i=l  Jtew(x)  - 

+  c^(S-x)  -  f  ♦.  ix^x-t)  tK  +  c'^(S-t)]  dt] 

+  L(S)  +c(S)  -c(S)  f  (|>,  ,,(x-t)  dt 
“  ■'teuCx)  '*'"2) 

-L(S)  f  <^,  ,,(x-t)  dt+  (n-2)  B(L(t))(S) 

~  •>^(x)  <«-!> -  - 

n**2  . 

+  (n-3)  B(c(t)){S)  -B(L(t))(S)  I  <>fi>(x-t)  dt 

-  -  “  "■  i=l  Jtea)(x) 

n-3  f 

-  B(c(t))  (S)  J  (x^)  dt  +  -{S,x) 

"  ■”  i^l^tewCx)  - n-2-- 

(n  >  2) 


(3.22) 


Next/  we  shall  study  the  asynptotic  bdiavioa:  of  fjj(x)  for  large  n.  For 

flO 

n  large  we  have,  because  of  the  convergence  of  J  .  (x-t)  to  \(»(x-t) , 

i*l 


y  (^  .  (x-t)  -►  4)(x-t) 
i«=l  (i) 

Moreover,  feom  Penerk  2.4,  P^,  i  ^  1,  is  a  cotpletely  continuous  linear 
operation  (\diich  might  vanish)  with 


iL  ^  If  2f  *  •  • 
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\*ere  r  and  e  are  positive  constants  independent  of  i.  Thus,  for 
large  n,  given  by  (3.20)  converges  to  a  function  D(S,x) 

vSiich  decreases  rapidly  as  n  becctnes  large.  Using  these  last  remarks 
in  (3.22)  we  get,  for  large  n, 

f„(x)  =  L(x)  +  f  t{>(x-t)  L(t)  dt  +  K 
“  ■'tewCx) 

+  c\[l  +  f  r|»(x-t)  dt]  +  c^(S-x)  +  L(S) 

^tew(x) 

+  c(S)  +  nB(L(t)+c(iy)  (S)  -  2B(L(t))  (S)  -  3B(c(t))  (S) 

-  B(L{t))(S)  f  t);(xHy  ^ 

""  ^tew(x) 

-B(c(t))(S)f  tt>(x-t)  dt  +  D(S,x) 

^tew{x) 

Or, 

fj^(x)  =  ng  +  u(x)  3hs:cr°  (3.23) 

v^e  g  =  B(L(t)4c(t))  (S^)  is  independent  of  n,  and  u(x)  is  a  function 

of  X. 

Now  fron  .(3.3)  for  large  n 

h  (x)  =  K  +  c''’(S-^)  +  f  (S) 
n  —  ~~  n  — 

«  K  +  c^ (S-x)  +  ng  +  u(S) 

xea 


ng  +  v(x) 


(3.24) 


I 


g  +  u(x)  =  L(x)  +  v(0)  tl  “  [  6(X“v)  djrl 


1 

.  1 

,  1 

f®®  r  T 

“  1  ^  +  1  '  C  jr  ,|,(3Qr)  ^ 

1  ,  ■'zewCx). 

1 

• 

1 

+  1  «(2)  6(x-v)  ^  xea° 

‘'Xew{x> 

(3.31) 

Let 

1 

1 

‘  ’  1 

H(x)  =  L(x)  -  [*"  (l><x-y)  ^  +  [  ■  j>(x-y)  ^  xea® 

'  ^0  '  'jjew(x)  1 

(3.32) 

■  • 

i(x)  =  [  ((.(x^r)  ^  , 

^Xew(x)  “  , 

(3.33) 

On  using 

(3.32)  and  (3.33)  in’  (3.31),  we  ctotain  after  transposing  g  to 

the  R.H.S 

•  ’  1 

■ 

1  '  :  ' 

u(x)  f=  -g  +  H(x)  +  v(0)  (1  -  $(x)) 

1 

i 

+  I  4)(x-y)  u(y)  ^ 

^Y.ei)(x) 

1  '  '  '  I  •  *  ; 

(3.34) 

Let 

1 

i 

A(x)  =  -g  +  H(x)  +  v(0)  [1  -  ^(jj)) 

(3.35) 

and  (3. '34)  beoones 

I 

■  f  ■  '  '  c 

u(x)  =  A{x)  +  oCx-y)  u{y)  ^  i  xeo  (3.36) 

—  “  Jtr _ M  — 


"  -I'cciV-n^ :7-;^ "y  : 
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Each  carponent  of  A(x)  as  given  by  (3.32)  and  (3.35)  is  bounded  over 
the  set  a)(S) .  Thtis,  A(x)  is  a  bounded  function  over  the  ocnipact  set 
w(S) . 

Prom  Theorem  2.1  the  solution  for  (3.36)  is  given 

u(x)  *=  A(x)  +  I  <l>(x-y)  A(y)  ^  ateo®  (3.37) 

“  ''2e(i»(x)  " 

«0 

vihere  \|»(t)  =  T  .  (t). 

“  n=l  “ 

Using  (3.35)  in  (3.37)  Me  get 

u(x)  »  -g  +  H(x)  +  v(0)  (1  -  i(x)l 

+  f  tl>(x«y)  t-g  +  H(^)  +  v(0)  [1  -  *(^))1  ^ 

^X?w(x) 

Thij  last  relation  is  true  for  all  xec®.  Set  x  =  S  and  we  get 

u(S)  «  -g  +  H(S)  +  v(0)  II  -  5(S)] 

+  [  il>(S-y)  I-g  +  H(^)  +  v(0)ll  -  5(j;)l)  ^  (3.38) 

''Xc«(S) 

Using  (3.28)  in  (3.38)  we  obtain 

v(0)  -  K  -  0*^8  =  -g  +  H(S)  +  v(0)  [1  -  ♦(S)l 

+  f  iIi(S-v)  I-g  +  H(v)  +  v(0)ll  -  ^  (3.39) 

T€(d(S)  " 

Let  us  denote  by 


-i'z^i,  ■-wz'^j^f\  -■.  .t^  »r/t; 
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*(n>(®'  ”  i  ♦/  ■  (S-Y)  d^ 


(n  >  1) 


(3.40) 


Note  ♦  .  (S)  =  ?  (S) 

(0  "  ” 


Roti  (3.39)  the  coefficient  of  v(0)  transposed  to  the  R.H.S.  is 


-i  +  (1  ~  $(S))  +  f  tl> (S~y)  (1  -  liy)]  ^ 
*'leu(S) 


-«S)  +  ?i,„,(s)  -  I  f  ♦(„,t£Z)  4®  &  0-«> 

m=l  n=l  '£ew(S)  '  ’ 


Fran  the  definition  of  5 .  (S)  (n  >  1) ,  the  g-oonvolution  pccpertiesr 

^(n)  ~  “ 


and  Renerk  2.3,  we  have  for  n  >  1 


^Zew(S)  'vcw(S)  'tc»i»(v) 


=  1  if 

Jteu(s)  heRj  (s-t) 


JtwCS) 

4>(n)  (S-t-y)  (j>(i[)  ^1  ^ 


f  4' is-y)  ^  =  I* ,  .  (S) 

Jloite)  '"♦i’  " 


(3.42) 


Using  (3.42)  in  (3.41)  wa  get 


-i@  +  1.  ♦(„)  ®  ♦(„«)  ®  =  ® 


n=l 


Simplifying  (3.39),  we  get 


g  U  +  f  \l»(S~y)  ^)  »  K  +  0*^8  +  H(S)  +  f  \l>(S-y)  H(y)  ^ 
^Y£(o(S)  ~  T®“lS) 


a 

1 


g  S3  M  {K  +  L(S)  +  [  t|)(S~y)  H^)  ^3 

^S~ye(D(S)  ‘  “ 


+  M  {  [  <|)  (S-y)  dy 

^ye(ij(S) 

+  [  t|>(S-y)  t  I  ((>(y-t)  t  dt]  dy 

^yew{S)  ^tew(y) 

-  f  c^y  <t>(S"y)  ^  } 


as  M  {K  +  L(S)  +  f  if<(S-y)  L{y)  ^}  +  c^y  -- 

^S-ycM(S) 

+  M  {  f  cV  <{)  (S~y)  dy 
^yea){S) 

+  I  f  ^ 

n=l  ^yew(S) 
n=l  ^yeuCS)  '  ’ 

Note  the  last  result  follows  from  the  definition  of  ,  the  convolution 
pcopeirties,  and  Lemma  2.1.  Upon  sinplifying  the  last  relation  and  sub¬ 
stituting  for  M,  (3.26)  follows. 

Q,  E.  D. 

As  given  by  (3.26)  the  stationary  expected  cost  per  period, 
g,  is  a  function  of  (0,S) ,  and  vd.ll  be  voritten  as  g(a,S) . 


■<  V  r^. 


‘?^r  T^-ytv*  ■«.»  ’»y4^  ■, 


OffiPIER  ly 

OHE  OFTlMIZmai  PROBLEM 

4.1  Introductaon 

In  this  chapter  we  shall  provide  necessary  and  svifficient 
conditions  for  the  exietenoe  of  an  optimal  (a,S)  policy  that  minimizes 
the  ei^ressicxi  for  g(a,S)  as  given  by  (3.26) .  The  unknowns  to  be 
determined  are  the  decision  variables  (i  =:  1,  2,  ...>  m)  and  the 
set  r. 

In  Section  4.2  xve  shall  first  determine  the  configuration  of 
the  set  r  and  then  reduce  the  minimization  of  g(a,S)  to  finding  optimal 
values  for  S  and  an  additional  variable  C,  vhere  C  will  be  identified. 

.  The  purpose  of  Section  4.3  will  be  to  determine  the  necessary 

and  sufficient  conditions  for  the  existence  of  the  pair  (C*,S)  that 
minimizes  g  (a  ,S) .  The  necessary  conditions  will  be  given  for  the 
gener^  m-oatinodity  inventory  problem,  vdiile  the  sufficient  conditions 
will  be  given  only  for  the  two-oarmodity  problen  with  seperable  L(x) . 

The  set  of  equations  satisfied  by  (C*,S^)  will  be  restated  in  terms  of 
a  real  valued  function  M(C*,x) ,  v^ere  M(C*,x)  will  be  defined. 

Section  4.4  will  give  a  geometrical  formulation  of  the  opti¬ 
mization  problem  in  terms  of  the  function  M(C*,x)  and  its  inherent 
properties. 

Finally  in  Section  4.5  the  linear  form  for  L(^) ,  the  conditional 
expected  holding  and  shox±age  cost  function,  will  be  considered. 
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4  >2  Characterisation  of  r 

In  this  sec±lon  vje  shall  first  (^varacterize  the  set  Tr 
r  ■  {2l|2S:n;  ZCS-x)  =  O),  by  determining  explicitly  the  relation  Z(S-x)  “  0 
iq>  to  a  constant,  nds  then  will  be  used  to  reduce  the  minimization 
problem  into  a  problem  in  differential  and  integral  calculiis. 

Theorem  4.1; 

Under  an  optimal  (a*,S^)  policy,  g(o*,S;*),  the  minimum 
stationary  total  expected  cost  per  period,  satisfies  the  following 
relation: 

g(o*,S*)  =  L(x*)  +  c\  x®er*  (4.1) 

vbere  ji  is  the  expectation  vector  of  the  random  variable  D  with  joint 
density  function  ^j^(t). 

Proof:  The  proof  is  based  on  the  definitions  of  h  (x) ,  xec,  and  f  (x) , 

'  ■  n  —  —  n  — 

jgEo  ,  as  given  in  Chapter  III.  At  optimality  we  must  have  for  x^eF* 
h„(x®)  =  f  (x«)  (n  >  1) 

We  note  that  this  last  relation  if  explicitly  established  will  determine 
the  configuration  of  r*. 

For  large  n,  ^  given  by  (3.1)  can  be  written 

f  (x)  «  L(x)  +  f  ((n-l)g  +  u(x-t)l  <^(1^  ^ 

^x-teh)(x) 

+  I  [(n-l)g  +  v(x-t))  (^(t)  ^ 


I 


-  I  [  (n~l)g  V  (x~t)  1  (t>  (t)  dt  ,  xeo*^ 

Vhere  g  is  the  stationary  total  expec±ed  cost  per  period  and  u(.)  and 

v(.)  are  well  defined  in  Chapter  III.  ■ 

•  1 

For  x®er*/  and  the  definition  of  u{x) ,  the  set  w(x*)  is  an 

I  t 

eirpty  set.  Hence#  from  (4.2)  we  can  write 

f^^®)  =  L(x®)  +  ((n-l)g  +  v(x®-t)]  4.{t)  x'eF* 

losing  (3-25)  and  the  asynptotic  esgaression  for  in  (4.3)  we 

I 

obtain 

ng  +  K  +  (S*-x°)  +  u(^)  ' 

«=  L(x®)  +  F  ((n-l)g  +  K  +  c'^(S*~x  -ft)  +  u(S*))  <l>(t)  ^ 
—  , 

I 

=  L(x®)  +  (n-l)g  +  K  +  c'^(S*-x®)  +  c\  +  u(Sf') 


siiiplifying  we  obtain  (4.1) . 


Q.  E.  H. 


The  result  of  Theorem  4.1#  as  we  shall  see,  later#  will  be 
of  great  inpcrtcince  in  finding  the  necessary  and  sufficient  ocaiditions' 
for  the  existence  of  eui  optiiral  policy  that  ‘ndnimizes  the  expressixm 
for  g(a#S)  as  given  by  (3.26).  An  explicit  expression  for  L(j<®)#  x®er*# 

I 

can  be  imnediately  derived  by  ocrtparing  (3.26)  and  (4.1)  to  yield 


L(x®) 


L(^)  +  f  L(S*-t)  ^(t)  dt 

^S*"tch)(S*) 

1+1  dt 

^S*-tew(S*) 


viiich  is  simil^  to  Iglehart's  result  [7]  for  the  single  cjomnodity 
case  operating  under  the  ( s»S)  policy.  For  notational  purposes  vje 
‘  shall  sometimes  represent' L(x“)  by  the  symbol  C*,  C*  refers  to  the 
value  of  the  total  expected  cost  per  period  excltiding  the  variable 
'  ordering  cost,  i 


feroof; 


T*  =  {x|xefi;  L(x)  -  c*  =  0} 


Prom  5?ection  2.2 

r*  =  {xjxen;  Z(S*-x)  =  0} 


Pron  Iheorem  4.1 


Hence, 


L(x)  -  C*  =  0 


I  • 

2(s*-x)  =  L(x)  -  C*  =  0  xer* 


and  (4.5)  follows. 


Q.  E.  D. 


let  us  make  the  transformaticm  x  =  S*~t;  then  (Chapter  II) 
the  image  of  the  set  o)(S*)  is  R(s*)  and  the  image  of  r*  is 
“  {i|t  >  0;  L(S*-^t)  -  C*  =  0  } 

The  sets  R(S*)  and  w(S*)  are  functions  of  C*  and  often  vdll  be  refers 

j  • 

to  as  R(S^,C*)  a«i  w(S*,C*). 


We  may  rewrite  the  esqjressicn  of  L(x") ,  x’ePr  as 
K  L(S*)  +  f  L(S*-±)  i|)(t)  dt 


The  relation  given  by  (4.1) ,  then,  nay  be  written 


L(S^)  +  I 


R(S*,C*) 


L(S*H:)  i)>{t)  dt 


Ftom  this  last  relation  it  nay  be  noted  that  to  mimnaze  g{OrS)  relative 
to  the  set  o  and  S  is  equivalent  to 

K  +  L{S)  +  f  L(S:;t)  i|>(t)  ^ 

Min  {  gj  (S,C)  = - : — 7 -  J 

i  +j  t|>(t)  dt 
-  ''r(S,C) 

vMch  is  a  problem  in  integral  and  differential  calculus.  Let 

M  =  (1  +  f  T|>(t)  ^  r*,  0  <  M  <  1, 

^R(S,C) 


then  the  minimization  problem  as  given  by  (4.8)  is  equivalent  to; 

min  {g,(M,C,S)  =  M  (K  +  L(S)  +  f  L(S-t)  i()(t)  dt]  } 

^R(S,C) 

S.T.  {gs(M,C,S)  =  M  U  +  f  tl>(t)  -  1  «  0  > 
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-A*  {1  -  M*  -  M*]  \[>(t)  dt  > 


rearranging  terms  we  obtain 


AG  =*  M* 


+ 


+ 


{L(^AS)  -  L(S*)  +  I  [L(S*-tAs-t)  -  A*]  \|r(t)  dt 

'  R(S*-4AS>C*4AC) 


I  [L(S*“t)  -  X*J  i|»  (t)  dt  +  AM  {k  +  L(S*4As)  -  A* 

^R(S^,C*)  - 

f  rL(S*4AS-t)  -  A*]  n»(t)  dt  } 

^R(S*4AS,C*4AC)  ” 


AA{1  -  (M*4AM)  -  (M*4AM)  rlf(t)  dt  }  (4.12) 

^R(S*4AS>C*4AC) 


If  we  denote  by  AR(S*,C*)  the  inoremental  set  between  R(S*^>vASfC*-tAc) 
and  R(S*,C*) ,  then  we  can  write 

f  [L(S*+As-t)  -  A^J  \(>(t)  dt 

^R(S*4AS,C*tAC) 


(L(S*~t)  -  A* I  i|»(t)  dt 

R(S^,C*) 


[L(S*+AS-t)  -  L(S*-t)]  t|»{t)  dt 
R(S*,C*)  “ 

+  f  rL(S*+AS-t)  -  A*)  rlt(t)  dt 

^AR(S^,C*) 


(4.13) 


and 
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J  CL(S*-«-/g-t)  -  X*]  t|;(t)  dt 

^R(S*+AS>C*4^)  - 

-  f  [L(S*->-AS“t)  -  X*]  t(»(t)  dt 

^R(S*,C*)  -  — 


+  f  [L(S*'t-AS~t)  -  X*]  dt 

ms*,c*)  - 


(4.14) 


Using  (4.13)  and  (4.14)  in  (4.12)  vte  obtain 

^G  =  M*  {L (5*4- AS)  -  L(S;*)  +  f  [L(S*+AS-t)  -  L(S*“t)]  i|»(t)  dt 

•'R(S*,C*)  -  “  — 

+  f  [L(S*4-/!S-t)  -  X*]  i|/(t)  dt  } 

J/SR(S*,C*) 

+  41  fK  +  •L'(S^+^)  -  X*[l  +  f  t(»(t)  dt) 

'R(S*,C*)  ”  "■ 

+  f  L(S*+AS-t)  ii>(t)  dt 

Jr(S^,C*)  “ 

+  f  .  [L(S*+AS-t)  -  X*]  i|;(t)  dt> 

^^(S*,C*)  ^  ~  — 


+  AX  a  -  (M*+ZM)  -  (M*+41)  f  4»(t)  dt  >  (4.15) 

^R(S*+ffi,C*+AC)  "" 


t®eiR(S^»C*)  and  let  A(^R)  denotes  the  area  of  the  incrennental  set 
/R(S*,C*) .  Then  from  Theorem  2.2 

f  tL(S*+^"t)  -  X*1  !(»(t)  dt ' 

•  ■'«(s*,c*) 

■  [L(S*4-^-t^  -X*)  \(»(t*)  A(m) 

t 

using  this  in  (4.15)  ws  obtain 


Let  t 
R*(-) 


A6  «  M*  {L(S*4AS)  -  L(^) 

» 

+  f  [L(S*+AS-t)  -'L(S*-t)l  i<>(t)  'dt 

^R(S^,C*)  ” 


+  [L(S*+AS-t.)  -  X*]  \l;(t®)  A(AR)}  +  AM  {K  +  L(S*+&S)  ‘ 


-  X*[l  +  f  t|»(t)  dtl  +'f  L(S*+AS-t)  il>(t)  dt 

•'r(s;^,c^)  •'R(ff»,c*) ,  “  ' 


+  [L(S*-«-AS-t®)  -  X*]  \|»(t")  A(AR)}  +  AX  {1  - 


(M^+AM)J 


^(t)  dt  } 


R(S*+^rC*+AC), 


jcr,*.  Expanding  L(f)  in  a  'mylcsr  series,  with  a  remainder, 

,  about'  {S*,C*);  the  expression  for  AG  as  given  in  (4.16)  becxmes 


m  3;i(S) 
AG  a  M*  {  J 


i«l  3S^ 


i  m  m.  3*L(S) 

X  " 


i=l,  j=l  3S^33j  ; 


AS^AS 
i  j 


I, 


m  3L(S-t) 

[  I 


AS. 


R(^,C*)  i*l  3S^.  &■§* 


,  m  m  3*L(S-t)' 

+  j  i  I  — - ASj^ASj 


1-1  j-l  3S,3Sh 


S*S* 


]  >(t)  ^ 


I 


r 


f 
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.  I 

I  .  •  » 

m  8L(S-t*) 

t  fL(S*-tJ  -  X*  +  I - AS. 

i-1  '  3S^  S«S*  ^ 


m  3L(S*-t) 

+  Y  -  At.}  t|>(to)  A(AR)>+  AM  {K  +  L(S*) 

i**!  3tj^  bsto  ~ 


m  3L(S) 

+  I  -  AS.  -  'X*  [1 

i=l  3Si 

j 

1 

1  '  1 

,  '  m  3L(S-tj  1 

+  j  [L(S*-t)  +  T  : - ini-  I  AS.l  i|»(t)  gt  ■ 

jR(S*,C*)  ,  i=l  3Si  ,S=S*  ^ 


'  m  3L(S-tJ 

+  lL(S*-tJ  -  X*  +  T - - 

,  -  m  3Si 


m  3L(S*-t)  '  ,  ' 

+  y -  ,  At.]  \;»(^)  A(AR))  +  .AX  U  -  {M*+AM) 

i«l  ,  3tji  t=to 


-  (M*+^)  f  ip(t)  ^  >  ' 

^R(S*+AS,C*+AC) 

» 

+  R*  (S* .G* .M* , X* S*+AS.C*-f AC rM*->-AMi'x*-fAX)  ’  (4.17) 


The  Abcfve  expression  for  AS  vdll  be  used  tc^  deteritdne  the  necessary 

\ 

and  sufficient  conditi^  for  '  (S,C)  to  have  a  proper  rela^ve 

1  »  ^  * 

mininun  at  (S^,C*),  i 


+  f  ^(t)  dtj; 

'R(S^,C*) 


1 
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Necessary  Conditions 

Wfe  shall  first  find  the  set  of  simultaneous  oonditicns 
that  (S*,C*)  vdll  satisfy  and  then  restate  these  conditions  in  terms 
of  the  function  M(C*,x)  vMch  will  be  defined  later.  Hie  following 
Zenma  is  needed. 

Lenina  4.1; 

At  optimality 

K  +  L(S*)  +  f  L(S*-t)  i|»(t)  ^ 

X*  - - -  C*  >  0 

1+1  ^ 

^R(Sf,C*) 


Proof;  The  partial  derivative  of  G{M,C,SrX) ,  as  given  by  (4.10)  r  with 
respect  to  M,  evalxiated  at  (M*,C*,S^,X*) »  yields 


3G 

3M 


(M*,C*,^,X*) 


«  K  +  L(S*)  +  L(S*-t)  \|»(t)  ^ 

^R(Sf,C*)  ^ 


-  X*  (1  +  f  «  0 

Jr(s^,c*) 

K  +  L(S*)  +  f  L(S*--t)  \)»(t)  ^ 

"  ''R(S^,C*) 


1  +  t(>(t)  ^ 

Jr  (§*,(:♦) 


ftan  this  last  result  and  (4.7) »  (4.18)  follows. 


(4.18) 


Q«  E«  D. 
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TBieorero  4,2: 

Ttie  optijtal  values  of  C  >  0  and  ^  C*  and  ^  /  vAiidh 
minimize  (S,C) ,  satisfy  the  following  set  of  equations: 

as,  s=s*  ^R(sf,c*)  as. 


C*  -  L(S*)  +  f  [C*  -  L(S*-t)3  tj>(t)  ^  «  K 

■'R(S*,C*) 


.  fin) 

(4.19) 

(4.20) 


Proof: 


The  partials  of  G(M,C,SA)  with  respect  to  (i  «  1,  2,  .  .  .,in) 
are  obtained  from  (4.17) 


aL(S) 


aS^  (M*,C*,S^,X*) 


=  M*(  — 

as. 

aL(s~t) 

i 

1 

dt 

aSi 

aA(/aR) 

+  [L(S*-to)  -  X*1  4i(to) 


J  «  0  (i  •  1,  2,.. ..a) 


where  the  point  tocFo*.  Since  by  (4.6)  and  (418) 

B  c*  «  L(S*-tJ  r  t^crp*f 

(4.19)  follows. 

TO  pcove  the  second  part  of  the  theorem  we  have  fecro  (4.18) 
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C*  [1  +  f  i,{t)  dtl  «  K  +  L(^)  +  f  L(S*-t)  ,p(t)  ^ 

^R(S*,C*)  ^  ■'R(S*,C*) 


And  this  last  relation  can  be  written  as 


C*  -  L(S*)  +  f  [C*  -  L(S*-t)l  \)>(t)  §t^K 

”  Jr(s*,c*) 


Q»  E«  D< 


Tbeorem  4.3: 

Let  the  function  M(C*,x)  satisfy  the  integral  equation 


M{C*,x)  =  C*  -  L(x)  + 

M(C*,x-t)  ((»(t)  dt 

(4.21)  1 

R(x,C*) 

i 

i 

i 

»  C*  -  L(x)  + 

4 

<|>(x-t)  M(C*,0  ^ 
u(x,C*)  ~ 

(4.22)  i 

Su{3pose  that  C*  and  exist/  vd^ch  minimize  (C/^)  subject  to 
C  >  0  and  S  >  0.  Tben,  if  >  0  and  C*  >  0,  (C*,S^)  is  a  solution 
to  the  set  of  equations 


M(C*,s;*)  *  K 


3M(C*,x) 

I  .  ■0  (i  •*  1»  2»  •  •  •/  w) 


Proof:  Define  the  functions 


.  3L(x)  f  3L(it"t) 

M.  (C*,x)  « - 11-  - \|i(t)  dt  (i  -  1,  2, 

^  "  3x^  'R(x,C*)  3X4^ 


(4,23)  i 

( 

! 

I 

«.2«)  I 


Hi  (4.251 


.  ft.  ^  ^  V  iL’sA 
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M(c*,x)  a  C*  -  L(x)  + 

f  rc*  -  L(x-t)]  t(;(t)  dt 

'r(x,C*)  “ 

(4.26) 

*  C*  -  L(x)  + 

f  [C*  -  L(t)]  T|>(x-t)  dt 

*^ci>(C*,x) 

(4.27) 

Let  X  =  S*  in  (4.25)  and  (4.26) 

then  ocmparing  the  results  with  (4.19) 

and  (4.20)  we  obtain 

/ 

M(C*,S*)  a  K 

(4.28) 

MjC*,S*)  =  0 

•  •  m  f  nO 

(4.29) 

Prom  Theoran  2.1  the  soluticai  of  (4.22)  is  given  hy  (4,27)  .  Htus 

M(C*,x)  a  M(C*,x) 

(4.30) 

Hence  fron  (4.28),  (4.23)  results  vMch  proves  the  first  part  of  the 
theorem. 

To  prove  the  second  part  of  the  theorem  we  can  write  from  (4.26) 

m(C*,x)  =  M(C*,x+Ax)  -  M(C*,x) 

\  «  C*  L(x4-Ax)  +  f  [C*  -  L(x4-Ax~t)]  i}>(t)  dt 

^R(x+Ax,C*)  “ 

-  IC*  -  L(x)  +  f  tC*  -  L(x-t))  4»(t)  dtl 

''r(x,C*) 

% 

If  we  denote  by  flR(x)  the  incronrental  set  between  R(x-fAx,C*)  and  R^,C*)* 
then  the  expression  for  ^M(C*,x)  can  be  written  as 


^(C*,X)  »  -  [L(xf  X)  -  LWl 


+  I  [C*  -  L(x+fl>c-t)]  ,j,(t)  dt 

JR(x,C*) 


+  I  [C*  -  L(x+^-t)]  a»(t)  dt 

JiP{x) 


f  CC*  -  L(x~t)]  dt 

Jr(x,C*) 


Let  t®eAR(x)  and  let  A(/S)R(x))  denotesthe  area  of  tf^(x).  -  Using  the 
Mean  Value  Theorem  for  multiple  integrals r  Theorem  2.2,  we  can  write 


/Sfl(C*,x)  =  -  tL(x+Ax:)  “  L(x)]  -  f  [L(x+AX-t)  -  L(x-t)3^(t)  dt 

iRCXrC*) 


+  [C*  -  L(xt^«)]  tl/CtO)  A(tfl(x))  (4.31) 


Ebqanding  L(x+^) ,  L(x-fAx-t”)  hi  a  Taylor  series,  with  a  remainder. 


R^  (xf Ax,x;C*) ,  about  x,  the  expression  for  tf1(C*,x)  becomes 


m  9L(x)  f  m  aL(x-t) 

^(C*,x)  =  -  T  - AX.  ”  I  - A?c.  4»(t)  dt 

i»l  3X£  ^  ^R(x,C*)  i*“l 


m  8L(x-t*) 

+  [C*  -  L(x-t»)  -  I  — ^(t«)  A(AR(x)) 

i-1  ax^ 


+  R*  (x»Ax»x>C*) 


As  to  -*«  0 ,  t ®  ^eFo  (x) .  Hence  lim  (L(x-t®)  -  C*)  -  0.  Dividing 

to+D 

^(C*,x)  by  toj  and  letting  tx.  -►  0  we  get 

3M(C*,x)  3L(x)  .  3L(x^) 

-  =  -  {  -  +  - 

3x£  3x^  iR(x,C*)  axj^ 

(i  «»  1,  2,  .  .  m)  /  (4.32) 


Ccxnparing  (4.25)  and  (4.32)  we  obtain 


3M(C*,x) 


=  -  Mj^(C*,x) 


Hence  frcm  (4.29),  (4.24)  follows. 


(i  ^  1,  2,  .  .  ni) 


Q.  E.  D, 


Sufficient  Conditions 

So  far  the  necessary  conditions  for  the  determination  of  the 
optinal  values  of  C  and  £  were  discussed.  To  study  the  sufficient  con¬ 
ditions  that  have  to  be  satisfied,  we  shall  assume  that  the  function 
L(x)  is  separable,  that  is, 

m 


L(x)  =  I  L  (x  ) 
i=l  i  i 


where  Lj>  (x^)  for  all  i  is  twice  differentiable.  Hgw,  since 
gs  (M*-f AM ,C*+£C , S*+AS)  ■  0,  and  since  the  necessary  conditions,  derived 
above,  require  that  each  derivative  of  6(M,C,S,X)  vanish  (4.17) 
reduces  to 


I 


I 


A<3  •  G(M*+AMrC*+AC,SJ^,A*+AA)  -  G(M*,C*,S*,A*y 


n  m  m  3L(S-t)  I  ,  ' 

"*'7.1,1,—; —  %is' 

j«l  i=l  3S.3S.  S=S*  ^  J 
1  3  ’ 

I 

■-  .  m  m  3L(S-t)  I  I  ,1 

■*  L/C*  .1  X  - -  '  ♦<«  at 

^R(S*,C*)  3=1  i«l  3S^3S.  ^S*  ^  j  — 


m  3L(S-tJ 
+  {L(S*-tn)  -  A*  +  I  —■ - ” 


i=*l  3S^ 


m  3L(S*-t) 

+  I  - =- 

i=*l  3t£ 


+  R*  (S^  ,C*  ,M* ,  A*  > S^+AS  ,C*+AC  ,M*+AM,  A*+AA) 

I 

vrfiere  t o  =  (t|g,  t2Q»  .  .  •,  t  )  cirg^r  RM«)  ’is  the  ratainder  terms 


(4.33) 


of  the  Taylor  series  ej^ansion,  A(AR)  is  the  area  of  AR(S^,C*),  and 

^'i^^i^'  I*”i(ti)f  (i  «  1,  2,  .  ,  m),  are  the  first  and  second  derivatives 

ofL^(t^). 

•  Since  L(S*-t^)  s»  A*  =  C*,  we  have  from  (4.33)  on  taking  the 
second  partial  derivatives,  for  i,  j  ■  1,  2,  .  .  m  » 


“7  J  -M*lV(S*)  +  f  L  "(S^.-t.)  ♦(t)  dt 


»(^)  I 


(4.34) 


I 


I 


3^G  ■ 


=  M* 


L  *  (S  \  tl'  ■  t  S  3A(AR) 


(C*,s*) 


(i?^) 


(4.35) 


3^ 

ls]lc  (M*,C*,S*,X*) 


3A{AR)  ' 

3C  {C*,S*) 


(4.36) 


Sc*"  (M*,C*,S^,\*) 


(4.37' 


Note  the  value  of  ^  is  not  the  sane  in  (4.34) ,  (4.35) ,  (4.36) ,  and 
(4.37) .  A  sufficient  ocndition  for  (4.9)  to  have  a  proper  minipiT"  a 
{M*,C*,S*)  is  that  tlie  matrix  A  given  below,  is  positive  definite. 
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1 

3*G 

3C3S, 

3»G 

3c3S2 

3*G 

3CSs 

m 

3*G 

3*G 

3*G  ... 

3*G 

acdSj 

&S,3S2 

3*G 

3*G 

3*G  .  .  . 

3»G 

3C9Sp 

3Sj  3S2 

3S7** 

3s.^s 
*  m 

• 

• 

• 

3*G 

3*G 

3*G  ... 

3*G 

dC3S 

m 

3S,3S 
*  m 

3S,3S 
^  m 

as^ 

m 

vrtiere  each  elenent  of  A  is  evaluated  at  the  point  . 


The  next  two  theorems  will  provide  the  sufficient  conditions 
for  the  twD-ocmnodity  problem  to  have  a  proper  minimiin  at  . 

Theorem  4.4; 

If  (C*,S*)  is  a  solution  to  the  set  of  equations  (4.19)  and 
(4.20),  then  (C*,^)  is  a  relative  minima  of  gi  (C,S)  if 

L^' (S£*“t^p)  <  0,  (i®l,2) ,  V^p  »  (tj prt2p)crp* 

L/(S.*)  +  f  L  "(S.*-t.)  i/>(t)  dt 

^  ^  ''R(^,C*)  ^  ^  ^ 

dA(AR) 


and 


+  L^*  (S^*-t^*)  \j>(t*) 


dSj 


(C*,S*) 


>  0 


(i-1,2) 


(4,38) 
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vi>ere  the  point  t*,  t*  «  (tj*,tj*)er,*,  has  dlffesrent  values  for  (i«lr2). 
ftpoofi  Fbr  the  two  ocrrtxJity  problem  we  mast  have 


3*G 

(a)  - >  0,  and  this  iinplies  firem  (4.37)  that 

3C* 


aA(AR) 

I 

ac 

•  I 

>  0 

(4.39) 

(C*,S*) 

a*G 

a*G 

I-gI  » 

1  -  — - 

- 1  >  0  and  from  (a)  tx.  xs  iitplies 

3c* 

as,* 

acasjJ 

a*G 

j 

1 

>  0 

(4.40) 

asj* 

Hence 

from  (4.34),  since  M*  >  0,  (4.38)  for  i  «  1  follows. 

3*G 

3*G  a 

*G  3*G 

(c)  {—  •  •—  -  ( - )*} 

as,*  ac*  asj*  acas, 


a*G 

a*G 

.  { _ . 

3*G 

3*G 

a*G 

3s,a', 

3c* 

3s,3s, 

3C3S, 

acas, 

o‘6 

3*G 

{ 

3*G 

a*G 

3*G 

•  - 

3c3s, 

acas. 

3Sj3s, 

3cas, 

as,* 

Upon  regiocping  terms  we  obtain 
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3*G  32G  a^G 

asj*  3C*  asi*  ■  ^  acasi 


3*G  a^G  a^G  a*G  . a^G 

+  2 - - - - - -  ,  -  jS 

3Si3Sj  acasj  acasj  ^2  as^as^ 


3*G  a*G  , 

- - ( - )*  >  0 

3S,2  gcaSj 


Upcjn  using  (a)  and  (b)  we  get 


a^G 

a^G 

3*G 

a^G 

( 

2  3*G 

3*8 

3*G 

3S,4 

30* 

aSi* 

acaSj 

38^382 

dCdS| 

3C3S2 

I 

Now  a  sufficient  condition  for 


(4.41) 


aA(AR)  I 

-  I  >  0 

3C  (C*,S*) 

is  that,  for  any  t,er,*,  L^'  (S^*-t^,)  <  0  (1=1,2) .  The  proof  proceeds 
as  follows: 


In  particular  if  A  =  (t„,tj,)cr,*.  Fig,  (4.1),  then 
C  =  I,,(S,*-t„)  t  I.,(S,*-t„) 

If  we  incteeient  C  by  je,  then  lor  B(tj^,tj^+Atjj)  and  D(t|j+4tjj,tj,) 


1 


I 
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(ii)  M(C*,x)  is  a  strictly  concave  function  of  x,  v^iere 
M(C*,x)  is  given  by  (4.21) 

Proof; 

Expanding  L(x+ax),  L(x-fAx-t°)  in  a  Taylor  series  vdth  a 
remainder  (x+ax»x;C*)  ,  about  x,  the  expression  for  in  as 
given  by  (4.31)  becomes 

2  2 

/SM(C*,x)  =  -  {  I  L.MxJAx.  +  I  Li"(x.)Ax.n  . 

i=l  1  X  1  1  1  1 


f  2 

-  H  Cl  L.'(x.-t.)Ax. 

^R(x,C*)  i=l  ^  ^  ^  ^ 


2 

+  I  L ."(x.-tJAx  *)  t|)(t)  dt)  +  {C*  -  L(x-t®) 

i^j_  1  •  1  1  1 - 


2 

-  <Xrti®)AXi}  t|<(t®)  A(AR(x)) 

+  R*  (x,x+^;C*)  (4.44) 

v4>ere  AR(x)  is  the  increnehtal  set  between  R(xtAxyC*)  and  R(x;C*), 

A(AR(x^))  is  the  area  of  AR(x) ,  and  t®EAR(x) . 

let  tj  =  (tio,t2o)ero(x).  Then  fron  (4.44)  we  can  write 

3*M(C*,x) 


{V(x.)  +  f  L."(x.-t.)  iKt)  dt 
^  iR(x,C*)  ^  ^  1 - 
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3A(AR(x)) 

+  L.'(x.-t.J  ,j,(t  )  - -  >  .(i=l,2) 

^  ^  3x.  (x,C*) 


(4.45) 


3*M(C*,x)  3A(AR(x)) 

- )  t|»(t  )  - ^  ijij  (i,j=l,2)  (4.46) 

3x^3xj  1  1  10  g 


Now,  M(C*,x)  is  a  strictly  concave  function  of  x  if  and  only  if  the 
Jfessienis  negative  definite.  Hence  we  must  have 


3*M(C*,x)  3*M(C*,x) 

(a)  - H.  <  0;  (b)  - 

3x  *  3x  * 


3*M(C*,x) 


3x,3x, 


3*M{C*,x)  3*M(C*,x) 


3x,  3x, 
1  2 


or,  fran  (4.45)  for  i=l. 


(a)  L /'(x.)  +  f  Lj"(x,-t,)  ii»(t)  at 

'R(x,C*) 


3A(AR(x)) 

+  L,'(Xi-tjo)  4'(t,)  -  >  0 

3xj  (x,C*) 


(4.47) 


3*M(C*,x)  3*M(C*,x)  3*M(C*,x)  a 

(b)  - =1-  .  - Z_  -  {  - }  >  0 

3x,®  3x2*  3xj3x2 


from  which 


r 


- Jf  I  k  ■  • 


3^M(C*,x)  ^  a2M(C*,x) 
3xi^  *  3X2"  "" 


Hence,  from  (a)  and  (4.45)  for  i  =  2,  we  have 


1-2  "(xi)  +  [  L2”(X2-t2)  li){t)  dt 

^R{x,C*)  - 


♦<£»>  (K,c»)  >0 
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(4.48) 


Set  X  =  s*  in  (4.47)  or  (4.48),  then  for  (i  =  1,2) 

Lj."(S^*)  +  f  V'^j*“t.)  \p(t)  dt 

''R(S*,C*)  -  — 


+  L.MS.*-t.„)  ,b(tD) 

1  1  lO'  TsT  (S*,C*) 


I 


>  0,  ^  =  (tio,  t2o)ero* 

Therefore,  condition  (4.38)  will  hold  if  M(C*,x)  is  a  strictly  ocncave 
function  of  x. 

Q.  E.  D. 

Uhfortunately ,  as  stated,  the  sufficient  conditions  are  not 
useable  in  practi<4^. 
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4.4  Gectnstric  Refopniilation  of  the  Problgn 

At  this  point  it  is  of  interest  to  use  the  inherent 
characteristics  of  the  function  M(C*,x)  in  giving  gecaietric 
reforraulaticn  of  the  problem  for  determining  the  local  minima  at  the 
point  (S*/C*)  in  E^. 

If  and  C*  satisfy  the  necessary  and  sufficient  oc  ?  ^icns 
for  a  proper  relative  minima,  then  as  we  have  seen  in  Section  4.5 
determining  and  C*  is  equivalent  to  constructing  the  function 
M(C*,x)  vhose  equation  is  given  by  (4.22)  with  the  followi:^ 
properties: 

(a)  for  x®er* 

M{C*,x®)  =  0 


3M(C*,x)  I 


(i  =  1.2) 


(b)  M{C*,x)  achieves  a  maximum  value  of  K  at  x  = 
i.e.,  M(C*,S^)  =  K 


Graphically,  this  reformulation  is  illustrated  in  Fig.  (4.2).  M(C*,x) 
is  plotted  for  various  values  of  the  paraneter  C*.  Ihe  function 


1 


:  a  =  a(K)  is  the  l6cus  of  all  points  S*  at  the  function  M(C*/X) 

attains  a  maximum  value  of  K  for  various  values  of  C*. 


4.5  Ihe  Function  L(x) 


Prom  now  on  we  shall  caansider  the  case  v^en  L(x)  has  the  form 


Kx)  =  J  L.  (X.) 

i=X  ^  ^ 

vAiere  each  L^(x^)  is  continuous  ai;^  twice  differentiable.  The  function 
L(x)  is  the  cost  charged  over  a  given  period  of  time  excluding  the 
ordering  cost  and  in  general  it  is  the  holding  and  shortage  costs. 
Consider  the  case  when  the  holding  and  shortage  costs  for 

I 

each  item  are  linear.  Let  for  item  (i  =  1,  2,  .  .  m) 


h.  =  Iftiit  holding  cost  at  the  end  of  the  period 

1  I 

I 

Pj^  =  Unit  shortage  cost  at  the  end  of  the  period 

) 

I 

Then  L(x)  is  the  total  e3?)ected  holding  and  shortage  costs  measured 
at  the  end  of  the  period  and  is  given  by 

m  ,  |Xi  p 

L(x)  ~  1  »».|'  (x.~t. )  (]) (t^ ^ .  ft- / . . • /t  )dt^ . . .dt. . •  .dt 

“  i=l  ^  JO  -'0  •'0  ^  ^  ^  i  « 


m 


ft^  f .  •  •  dtj  • .  .dt^ « . 


(4.49) 


I 
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If  and  $^(*)  denote  respectively  the  itargiiial  density  function 
and  marginal  distribution  function  of  the  demand  for  item 
i  (i  =  1,  2^  .  .  m)  in  a  given  period,  then  from  (4.49)  we  obtain 


m 

L(x)  =  I  L.  (x.) 
i=l  ^  ^ 

L^(Xj^)  denotes  the  marginal  e:q)ected  holding  and  shortage  costs  for 
iten  i  measured  at  the  end  of  a  period.  Fran  (4.51)  it  is  clear  that 
the  function  is  twice  differentiable  at  all  points  for  vMch 

(x^)  is  continuous.  In  fact  for  (i  =  1,  2,  .  .  .,in) 
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Li’(Xi)  =  (hiHp£)4»£(Xi)  -  Pi 
“  (hi+Pi)4'i(Xi) 

and  L^Cx^)  is  stxictly  cxsnvex  for  all  points  for  v»hich  ((i^Cx^)  y  0. 
Henccr  L(x) ,  being  the  sum  of  strictly  convex  functions,  is  itself  a 


strictly  convex  function. 


C3fflPTER  V 


OClMPUTATIC»iaL  ASPECTS  OP  TOE  PROBLEM  -  AN  EXAMPLE 

5.1  Introducbion 

In  this  dicpter  the  analysis  shall  be  restricted,  to  the 
oonputational  aspects  of  the  optimization  problem,  for  the  special 
case  of  a  tv^o-oonmodity  problem  v^ere  the  demand  for  the  items  dber/s 
the  e^^xjnential  distribution,  and  the  holding  and  shortage  costs  are 
linear. 

In  Section  5.2  we  will  present  the  two-ccmodity  problon 
xmder  consideration  and  obtain  e3?)licit  e^^ression  for  the  function 
LfCx^) ,  (i  =  1,2) ,  as  given  by  (4.51),  and  its  first  and  second  derivatives. 

Ihe  purpose  of  Section  5.3  will  be  to  determine  analytical 
egressions  for  the  set  of  equations  used  to  determine  the  c^timal 
policy  parameters.  The  integral  equation  (4.22)  will  be  converted 
into  a  partial  differential  equation  of  the  second  order.  Riemann 
method  will,  be  used  to  solve  this  boundary  value  problen  at  the 
point  S . 

Finally  in  Section  5.4  a  numerical  exanple  will  be  considered. 
Nunerical  methods  will  be  used  to  solve  for  S^*,  S2*,  and  C*. 
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,1 


5,2  Exanple 

Consider  a  tvro-<3csmodity  inventory  control  problon  v^iere  the 
demand,  {Dj},  for  the  items  over  a  sequtoioe  of  periods  (j  -  1,2,...), 
is  assxxned  to  be  ind^iendently  and  identically  distributed  oontinijous 
randan  variables  with  joint  density  function 


=  4>(t)  =  XiX2e"^i^i"^2t2 


Ihe  marginal  probability  distributicHi  function  is  given  by 


♦i(t^)  =  <|)i(t)  dt  =  1  - 


Fran  (4.51)  we  can  write 


L^(x^)  =  (h^  +  p^)  I  ^  ‘^^'i  ^i^^i  ”  ^i^ 


=  +  Pi)  I  ^  (1  -  '3ti  +  Pife  -  *i) 


=  (h,  ^  p.)  (x.  ^  -  ji.)  +  Pi(i  -  xj 


(h.  +  p.)  ,  h. 

=  L^.  .-..y,  +  h.x.  - 


1  i  X. 


(i  =  1,2) 


(5.1) 


On  taking  the  first  and  second  derivatives  of  (5.1) ,  we  get 


L.'  (x.)  =  -(h.  + 


P.) 


e'^i^i  + 


(i  =  lp2) 


(5.2: 


J 


“XiXi  : 

l5^?'(x  )  =  X£(h^4p^)  e  >0  (i=l,2) 


Using  (5.1)  in  (4.5)  we  get 


r  I  L  ^1*^1  “XiXi  h^+p2  -X»x, 

r*  -  ixjx^En;  hjX,  +  hjXj  + - e  + _ e  ^  * 

Xi  Xi 


‘  105 


(5.3) 


ha 

=  C*+  —  + _ } 

X»  Xa 


(5.4) 


Prcm  Theorem  4.2,  Sj*,  s^*,  and  C*  are  the  solutions  for 
the  set  of  equations  given  by  (4.19)  and  (4.20).  Using  (5.1)  and 
(5.2)  in  (4,19)  eind  (4.20)  we  get 

1 

hi  ~  t^i+Pi)  ®  ^  ^  +  f  [hi  -  (hi-fpi). 


-Ai(Si»-ti) 

e  J  ,j,(t)  ^  =  0 


(i=lr2) 


2  (hi+Pi)  -AiSi*  ^  hi  ' 

C*  -  J  [ - e  +  hiSi*  -  —  ]  +  f  .  ic* 

H  '  'R(s*.-c*) 


2  (hi4pi)  -Ai(Si*-ti) 

’  ii  ~7~  ® 


(5.5) 


—  1  »(>  (t)  dt  =  K 
X; 


(5.6) 


I 


1 


V'(t)  =  I  ^  (t) 

n=l 


neon  (5. 5)  we  can  write 


f  “^i^i 

I  (h4+Pi)e  ih(t)  dt 

.  -XiSi*  ^R(S*,C*) 

[^i  “  e  i  1  =, - - “7 - i— ^ 

,  1  +  .  ’I>(t)  dt 

■'R(S*,C^) 


f  >1 

e  t|>(t)  dt  I  : 

,  1  ^R(S*,C*)  .  hi  . 

Si*  =  —  In  [  r—. - ; - : - - r-  1  (i«l,2) 


—  jji  I  f  I 

Xi  .  ^{t)  dt 

•'r(S*,C*)  . 


hi+pi 


The  cenplexity  of  finding,  an  analytical  expression  for  ip  (t)  dictates 
abandoning  the  direct  solution  of  (5.5)  and  (5.6)  for  Si*,S2*»  and  C*. 
Now  frem  l^ieorem  4.3  the  values  of  S*  and  C*  are  the  real 

i 

positive  soluU.on':  to  the  set  of  simultaneous  equations  givm  by  (4.23) 
and  (4.24)*  i.e.  ' 


9M(C*,x) 

M(C*,S*)  =  K;  ^ -  =  0 

9x .  x=S*  ' 


(i=lr2) 


vdiere  M(C**x)  is  given  ^  (4.21) 


nv*.?" 


:T* 
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To  solve  the  set  of  equations  given  by  (5.7) ,  we  shall  first 
find  an  expression  for  M(C*,x)  and  then  use  this  expression  to  deter¬ 
mine  Sj*,  Sj*,  and  C*.  One  way  to  solve  for  M{C*,x)  is  to  cxanvert  the 
integral  equaticm  (4.22)  into  a  partial  differential  equaticai  of  the 
second  order. 

Iheoreni  5.1; 

In  E*  if 

-Xjtj  -  Ajt, 

(|>(t)  =  e  (X^  >  0,  tj^  >  0;  i=l,2) 

and  L(t)  =  Lj(tj)  +  L2(t2),  where  i=lf2,  is  given  by  (5.1), 

then  M(C*,x)  satisfying  the  integral  equation  (4.22),  i.e. 

M(C*,x  )  =  C*  -  L(x)  +  f  <|)(x-t)  M(C*,t)  dt  (5.8) 

^a)(x,C*) 


is  the  solution  for  the  boundcury  value  prdblem 


3^M(C*,x)  3M(C*,x)  9M(C*,x) 

I  I  ^  ^2  I  ^  ^  ~  ^*^2  ^ 

SXjSXj  3Xj  axj 


and  for  x®er* 


M(C*,x®)  =  0 


(5.9) 


3M(C*,x)  I 
3Xi  3C^® 


-V'v 


(h^-«P^)e 


-X,x,® 


-  h. 


(i=l,2) 
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Proof; 

B«?fore  deriving  the  partial  differential  equation,  we  will 
iiota  that  for  i,j=l,2 


3<l>{t) 

— ^ — =  -Xi((>(t) 

ati 


-  =  ).iX4(})(t) 

at^atj  ■' 


The  partial  differential  equation  will  be  derived  directly  from 
the  integral  equation  using  again  the  method  of  variable  increment  on 
M  as  in  Chapter  IV. 


AM  =  M(C*,xtAx)  -M(C*,x) 

=  C*  ~  L(x4-Ax)  +  f  (x-f Ax-t)  M{t)  dt 

*a)(xtAx,C*) 


-  [C*  -  L(x)  +  [  (|)(x-t)  M(t)  dt] 

'(dCx^C*) 


adding  and  subtracting  4>  (xtAx-t)  M(t)  ^  to  the  right  hand 

^(o(x,C*) 


side  we  get  after  simplification 


AM  =  -  {L(x+Ax)  ^  L{x)}  +  f  [4' (x+Ax-t)  -  4»(x-t)]  M(t)  dt 

~  ^o)(x,C*) 


(5.10) 


(5.11) 
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‘-J’ 


no 


,  2  2  3»(i)  (x-t) 

+  y  Y  y  - Ax*AXi  1  M(t)  dt  +  {({»(x-t®) 

j=i  i=i  ax^axj  ^  ^ 


2  3(i)(x-t°) 

+  y  - AX.}  M(t®)  A(Aa)(x))i  -f  (xfAx,x;C*)  (5.14) 

i=l  3x^  **  “  “* 


v^iere  (x+AX/k;C*)  deriotes  the  remainder  term  of  the  es^Jansion.  On 


using  (5.2),  (5.3),  (5.10),  and  (5.11)  in  (5.14)  we  get 


“XiXi  -XiXj 

AM={(hi+Pi)e  -hi}Axi-{(h2+P2)e  -  h2)  Axa 


“XjXj  “^2^1 

•-Xj(hj+Pj)e  Axj*  -  Xj{h2+p2)e  Ax2 


+  J  (  “X  i(i(x-t)Ax.  “  X  (|»(x“t)Ax2 
fc)(x,C*)  . 


+  X,^4i(x-t)Ax*  +  X^X2(t>(x-t)Ax^Ax,  +  X,^(^(x~t)Ax,^l  • 


^  ^‘^(x~t°)  -  X,(t»(x“t^)Ax,  -  X,(^(x--t°)Ax,}  • 


M(t^)  A(Au(x))  -f  (x-fAx,x;C*) 


Figure  (5.1)  shews  geonetrically  the  incxanental  areas  due  to  positive 
increments  in  Xj  and  Xj. 

New  fron  (5.15)  we  have  for  i=l,2 


.  3M(C*,x)  -X.x.  f 

- ^  =  (h.+p,)e  -  h.  -  X.  J  (J>(x-t)  M(t)  dt 

3Xj^  ^  ^  ^  ^  '(i)(x,C*) 


(5.15) 


3A(A£i)(x)) 
+  M(^) - z- 

3}^ 


^eAu(A3^) 


or  eaqplicitly 


3M{CSx)  -X,x, 


3x, 


A|A.|  0 

=  (h,+p,)e  -  h,  -  X,  J 


w(x,C*) 


<^(x-t)  M(t)  dt 


+  <^(0,Xj-t2o)  M(Xj,t2,) 


3A(Aa)(x)) 


where,  from  Pig.  (5.I),  t^jE  [Xj^Xj]. 


3M(c*,x)  -XjXj  , 

-  (h2+P2)e  -  hj  “  X2j  (t> (x-t)  M(t)  dt 

‘'w(x,C*) 


3x, 


+  <I'(xrt,,,0)  M(t,,,X2) 


3A(Aw(x)) 

Ax, 


where,  f ran  Pig.  (5.1),  t,oe[x,,x,I. 


3*M(C*,x) 
3x  j3x2 


=  X,X,f 

1  2j 


a)(x,C*) 


(^(x-t)  M(t)  dt 


3A(Aa)(x)) 

Xj((.(0,X£-t2e)  M(Xj,t2o)  - ^ 


3A((!6)(x)  ) 

-  X'jV(X,-tjj,0)  M{t,5,X2)  - =- 

3X2 


(5.16) 


(5.17) 


Xj+Axj 
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+  J  [c*  -  h|X,  -  hjXal  . 

w(C*,S) 


-X,(Sx-X|)  “  Xj(Sj-X2)  _ 

e  I^[2»'XjX2(Si-Xj)  (Sj-Xj)]  dxjdxj 


v^erc  Ijt  ]  is  tlie  modified  Bessel  function  of  the  first  kind  of  zero 
order. 


Proof: 

3h  solving  the  boundary  problem,  which  is  a  linear  hyperbolic 
partial  differential  equation  of  the  second  order,  we  will  use  a  method 
due  to  Ri-emann  (see  P9 ) .  The  value  of  the  function  at  tlie  point  S  is 
given  by 


1  (  3m(C*,x) 

K(C*,S)  - - {  U(S,x)  - TLdx, 


+  X|X2[c*  -  hjx^  -  h2Xj'J  U(S,x)  ^  ) 

^w(C*,S) 


v^iere  U(£-x}  is  the  -*een's  function  for  the  problem  satisfying  the 
following  boundary  problem 


3^U  9U  3u 

- -  X. - -  xi - =  0 

SxjSxj  9x,  9x2 

9U 

- *=  X.U  vten  X-  =  S, 

^  »  2  2 


(5.19) 


(5.20) 


(5.21) 
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8U 

3X2 


=  A2U  viien  Xj  =  S| 


U  =  A1A2  when  Xj  =  Sj,  X2  =  S2 


(5.21) 


A  solution  for  (5.21)  is  given  by  (see  [25]) 


U(x,S)  =  A1A2  er^l<Sj-Xi)-A2(S2-X2)  lo  X2)  ] 

•M*  •>»  ^  ^ 

(5.22) 

using  in  (5-20)  e^^jression  (5.22)  and  the  valiie  of 
^HC*,x) 


as  given  in  (5.9),  \ie  obtain  (5.19) 


Q.  £.  D. 

5.4  Numerical  Exanple 

In  this  n\jn2rical  exanple,  we  shall  make  use  of  the  parameter 
values  considered  by  Sivazlian  [13] .  Let  us  consider  the  case  vben 

Ai  =  A2  =  1 

Pi  =  P2  =  20 

hi  =  h2  =  1 


■i^ 


K=5.0 

Using  these  values  in  (5.1),  (5.2),  (5.4),  and  (5.19)  we  get 


L  (x  )  =  +21e  ^  +  X£  -  1  (i=l,2) 

X  i 

Li'(xi)=-21e  ^  +  1  (i=l/2) 


-Xi  “X, 

r*  =  {2L|^i  +  21e  +  Xj  +  21e  "  =  C*  +  2  } 

and  at  the  point  S  =  (8^,82) 


,  ■"(Sj“Xj)  ”* 

M(C*,S)  =  e  Ij2v/(S^-x^r(8^~xp3  . 

J  p* 


{21e  ^  -  1}  dx,  +  [  (C*-x,-x  )  . 

^u)(C*,S) 

~(S.-Xj)  -  (8  -x,)  _ 

e  Ip[2/'.S,-Xj)  (Sj-Xj))  ^ 


Fig.  (5.2)  illustrates  geonetrically  the  configuration  of  the  set 
a)(C^S)  for  the  case  when  C*=  8.  To  determine  the  values  of  8^*,  8^*, 
and  C*,  we  have  to  solve  the  set  of  simultaneous  equations  given  by 


(5.7),  v4iere  M(C*,x)  | 

x=8 


is  given  by  (5.23). 


The  conplexity  of 


obtaining  an  analytical  expression  for  M(C*,S^)  is  evident  and  therefore 
numerical  results  are  s?*xght. 


« Jii  '  '"  -  '  >  I  h  I  iH  I  ji.  ■TirrrrprTrn’TT^  .1  1  ii  1  '■i*'"*  ■■ 
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Nurtierical  Analysis; 

In  carrying  out  the  numerical  analysis  for  this  particular 
exanple,  the  followirig  points  were  taken  into  consideration: 

(i)  at  optimality  Si*  =  S2*  , 

(ii)  for  a  given  value  of  C  ,  T  is  syimetric  about  the 
line  Xi  =  X2  ‘  , 

I 

(iii)  the  results  obtained  by  Sivazlian  [13] ,  table  1,  for 

i 

this  particular  exanple  under  different  (a,S)  policies 

(iv)  Sivazlian  [13]  shows  that  . 


[  [  io[2/uiu^  dui  du2 

Jo  '0 


=  a{l  -  [e”^®  Io(2a)  +  e'’^^i(2a)  ]}  (a  >  0) 

i 

(v)  Sivazlian  [13]  further  sinews  that 


f  f  u.  e“  Iq  [2/ui U2 J  duj  du2 

Jo  Jo  ^ 


=  {1  -  Ip  (2a)  +  e“^®  Ii(2a)]} 

■  (i  =  1,2)= 


(vi)  for  this  particular  exariple,  M(C*,x),  as  giv^  in  (5.8) 


is  a  non-decreasing  function  of  C  . 


I 


I 


Taking  the  above  points  into  consideration/  we  initiated  the 

i  t 

search  procedure  by  choosing  C*  =  8.2  ^d  8^=82=  3.8.'  In  confuting 
the  double  integral  in  (5.23)  for  a  particular  value  of  C*,  S,  ard  5,, 
the  region  («)(C*,S^)  is  partitioned  into  four  smaller  jfegions  as  shown  in 

I  • 

Tig'.  (5.3).  Th''  crv^'^dimtes  of  the  point  Pj  aro  (Pj-3  .PfR^-l.S) . 

'•I  ■  ,  .  ’ 

'  (1)  ■  Using  (iv)  and  (v)  we  have 


f  -(S,-x,)  7  (SjTXj) 

(C*-x,-X2)e  Ij 


-Xj)]  ^ 


-3.6  -3.6 

(C*-Si-82)  [1.8{1  -  [e  Ip(3.6)+e  Ii(3.6))}] 


1.8[1  +  3.6]  .  -r3.6 

+  2 - {1  -  (e  Io(3.6) 

.  2  >  ' 


-3.6  (1.8)2 

+  e  I,  (3.6))) - 

,  2 


=  1.07109288  (C*-Sj-S2)  +  1/68702725 
-3.6  -3.6 

‘  v^ere  the  values  of  e  Iq(3.6)  and  e  Ii(3.6)  were 

I 

taken  frctn  Bessel  function  tables  [19] . 

(2)  The  double  integral  of  the  function  , 


— (8^— Xj)  —  (S2—X2) 

(C*-Xi-x^ e  I,  r2/(Si-Xi)  (S2-X2)  ]  / 


1 
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vAiich  is  defined  over  the  closed  subregions  A2»  A3 
and  At,  is  formulated  hy  cutting  each  subregion  into 
very  small  squares  each  of  area  (i  =  2,3,4) , 
vAiere  any  two  of  these  squares  have  at  most  one  side 
in  oornnon.  Hence  for  (i  =  2,3,4) 


(C*  -  xi  -  X2) 


(Si-Xi)  (S2-X2)]  ^ 


=  liin  Y  (C*  -  xi  -  X2  )  >^®2-X2  ) 

Ahi-’-O  j=l  j  j  J  J 


Io[2/(Si  -  Xi.)  (S2  ~  X2.)  ]  AA. 

3  3  ^ 

Valero  (xi .  ,X2 ,)  is  any  point  in  A.  and  n.  is  the  maximum 
3  3  13. 

number  of  squares,  of  area  AAj|^  in  Aj^.  synnetry,  the 
integral  over  A4  is  equal  to  the  integral  over  A2. 

(3)  The  line  integral  of  the  function 


^j-(Si-Xi)-(S2-X2)j^j2/(Si  -  xi)  (S2  -  X2))  •  {21e~^  -  1} 

which  is  defined  over  the  smooth  arc  F*  is  formulated 
by  dividing  the  arc  F*  into  N  arcs  hy  inserting  (N  -  1) 
points  Qj|^{xij:,X2i) .  Thus, 


g-(Si-Xi)-(S2-X2)j^j2  (Si  -  xi)  (S2  -X2)  ){21e“^l-l]}  dxj 
F* 
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=  lim 

Axii^O 

Ax2i-^ 


I  e“ - (S2~X2i) j2/(Si-Xii) (S2-X2i) ] {21e“^li  -  1}  Axji 
i=l  • 


Nuirerical  Results; 

To  obtain  numerical  results  /  a  ccnpatational  algorithm  in 
Fortran  IV  (appendix  5. A)  was  written  tc  solve  for  the  optimal  policy 
variables.  The  algorithm  proceeds  in  the  follcwing  steps: 

Step  0  -  Choose  initial  values  for  C  and  tlie  incraiient  in  C  ,AC  ; 

-  Choose  initial  values  for  S /  S  =  So ,  and  an  increment  in  S , 
AS.  The  range  of  values  of  S,  Sq  £  S  £  S®,  -for  a  given  value 
of  C*,  are  chosen  to  assure  that  M(C*,S)  achieves  a  maxirauin 
in  that  range; 

-  Qioose  AK  as  an  upper  limit  on  K  *•  M(C*,S) . 

Step  1  -  Ccmpute  M(C*,S) . 

Step  2  -  Check  the  difference  between  M{C*,S)  and  K,  then 
a.  If  the  difference  is  greater  than  AK,  set 

AC*  =  ^  ,  C*  =  C*  -  AC*,  S  =  So  and  go  to  Step  1. 
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b.  If  the  difference  in  absolute  value  is  greater  than  M., 
go  to  Step  3. 

c.  If  tlus  is  the  first  time  the  difference  in  absolute  value 
is  less  than  AK;  set  K*  =  M(C*,S)  /  S*  =  S/  Cj  =  C*  and  go 
to  Step  3.  Otherwise,  go  to  Step  3. 

Step  3  -  Set  S  =  S  +  If  the  new  value  of  S  is  in  the  range  of  values 
of  S,  go  back  to  Step  1.  Otlierwise,  if  K  -  M(C*,S)  >  AK  over 
all  the  range  of  values  of  S,  set 

Ac* 

AC*  =  __  ,  C*  =  C*  +  AC*,  S  =  S 0  and  go  back  to  Step  1;  else 

the  algoritlTO  converges  on  the  optimal  values  of  C*  and  S*. 

For  tlie  particular  exanple  we  are  considering,  when  K  =  5.0,  we 
initiated  the  iterative  search  procedure  as  follows:  Tlie  values  of  C* 
were  chosen  over  the  range  of  values  [8.22,  8.27]  in  increment  of  0.01; 

Si  =  S2  Values  ranged  over  the  closed  interval  [3.8,  4.08]  in  increment 
of  0.02.  The  number  of  "steps"  used  in  the  numerical  integration  Nj 
were  selected  to  be  60  and  100.  For  the  case  Nj  =  100  we  have 

C*  =  8.2625  and  S*  =  S«  =  4.02 

M(C*,  S*)  =  5.00C94  and 


9M(C*,S^) 


S 


S* 


0  for  (i 


1,2) 


wiiich  ca^paro  favorably  with  the  results  obtoined  hy  Sivazlian  [13] 
vdien  foliating  different  optimal  policies. 
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Graphically,  the  cxxiputed  values  (Appendix  5. A)  of  the  function 
M{C*,S)  cire  plotted  in  Figs.  (5.4)  and  (5.5)  for  the  cases  when  =  60 
and  100  respectively.  As  an  experimental  observation,  it  my  be  noted 
that  as  the  number  of  "steps"  used  in  carrying  out  the  numerical 
integration  increases  the  curves  become  smoother. 

'  Ihe  aocuracy  of  the  oenputed  results  is  subject  to  two  sources 

of  error,  truncation  and  roundoff.  In  cenputing  the  modified  Bessel 
function  Iq  (x) ,  only  the  first  ten  terms  of  the  convergent  series  were 
sunrned.  However  it  was  found  ejperinentally  that  the  roundoff  error  is 
the  factor  of  significance  and  it  is  not  ea;^  to  examine  it 
analytically,  since  it  is  a  function  of  the  oonputcr  used  and  the  size 
of  the  "steps"  used  in  carrying  out  the  numerical  integration. 


,.  .  .  '  CHAPTER  VI 

CmCLUSIONS  AND  REOCSV5MENDATICNS 

1  '  I 

6i.l  Ctotnclusions  ' 

I  • 

A  Stationary  multi-cqnnociLty  inventory  problem  with  periodic 

i 

review  has  been  formulated  and  studied  when  a  dyadic  {a,S)  policy  is 
used.  At  tlie  beginning  of  each  of  :a  sequent  of  periods  of  time  the 
stock  level  of  eac^  item  is  reviewed  and.  a  decision  to  order  or  not 
to' order  is  made.  Ihe  cost  elsnents  that  affected  the  ordering  deci- 

'  I 

in  each  period  are:  the  ordering  costs,  the  holdiiig  cost  and  tJie 

.shortage  cost.  Tne  ordering  ajsts  are  assumed  to  be  cenposed  of  a 
1  ' 

purchasing  cost  which  is  proportional  to  tlie- quantity  of  each  itan 

I 

1 

ordered  and  a  single  set  up  cost  which  is  independent  of  the  quantity 
ordeiW.  The  holding  and  the  shortage  costs  are  assumed  to  be  diarged 
on  the  basis  of  the  stock  leve^  at  the  end  of  the  period.  Demand,  for 
the  items,  in.eacli  period  of  time  is  described  by  a  continuous  random 

I  ’  ~ 

vector,  with  a  joint  density  function,  independently  distributed  frem 
period  td  period.  '  Imnediate  delivery  of  orders  and  conplete  backlogging 
of  all  unfilled  denumds  ctre  assumed.  ■ 

Tlie  ihventory  systan  just  doscrited  was  treated  by  many  researchers 
under  different  ootimai  (a,S)  policies.  This  research,  as  cettpared  with 
wha(t  has  been  done  .in  tliis  area,  .is  more  general'  in'  •the  sense  no  cissunp- 
tions  were  made  about  the  configuration  of  tlie  ordering  region  and  no 

I  '  ' 

specif ici  joint  density 'functions  were  considered. 


! 
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The  m-dimsnsionnl  convolution  operation  wliich  was  introduced 
in  Chapter  II  and  used  in  the  study  is  a  generalized  concept  of  the 
ordinary  convolution  operation.  For  the  case  when  I’,  the  hiTpersurface 
that  si±idividGS  the  ordering  region  and  the  not-ordering  region/  is 
adinissible,  the  g-convolution  of  any  function  reduces  to  the  ordinary 
convolution  iuid  all  the  g-convolution  operation  properties  with  the 
Cenmutative  hold.  The  g-convolution  properties,  in  pcurticulcur, 
were  used  in  tlie  solution  of  an  m-diiincnsional  integral  equation  of 
the  renewal  type,  which  is  similar  in  form  to  the  equation  solved  by 
other  rcsearcliers  in  tliis  area  but  with  no  restrictions  on  the  con¬ 
figuration  of  r. 

The  analysis  used  in  deriving  tlie  analytical  expression  for  tl^e 
stationary  level  ei'pected  cost  per  period  is  similar  to  Sivazlian’s  [113 
approach  for  the  one  ooiTmodity  probleiT)  when  operating  under  £;  station¬ 
ary  policy  of  tlie  (s,S)  type.  The  asi^iptotic  results  for  the  prc±)lan 
were  deduced  by  appealing  to  the  results  of  Yosida  and  Kakutan  [22]  in 
the  theory  of  linear  operations. 

At  optimality  it  was  shown  that  the  set  T*  that  subdivides  tlie 
policy  space,  fi,  is  given  by 

r*  =  {x|x  e  Q',  L(x)  -  C*  =  0} 

v^herc  C*  is  the  minijnum  value  of  the  stationary  total  cxj[xjcted  cost 
per  period  excluding  tiie  variable  cost,  and  L(x)  is  the  conditional 


expected  holding  and  shortage  cost  function.  By  using  this  result, 
the  minimization  problem  Wcis  reduced  to  finding  the  decision  variables 
Sji  (i  =  l,2,...,m)  and  C  tliat  minimize  the  stationary  total  expected 
cost  per  period.  S*  and  C*  were  given  as  the  real  positive  solutions 
to  the  set  of  simultaneous  equations 

M(C*,S)  =  K  (6.1) 

3M(C*,x) 

— r- -  *  =  0  (i  =  l,2,...,m)  (6.2) 

v^iidi  are  similar  in  form  to  the  set  of  equations  used  by  Sivaz^'^an  [11] 
in  determining  the  optimal  values  s*  and  Q*  =  S*  -  s*  for  the  one 
oomrodity  prcblon.  M(C*,x)  satisfies  the  integral  equation 

/  M(C*,x)  =  C*  -  L(x)  +  I  M(C*,x-t)  (J)(t)  dt  (6.3) 

•'R(x.C*) 

in  viiich  ^  is  the  joint  density  ftinction  of  the  dsnand.  For  the  case 
of  a  two-ocranodity  problem  \\here  the  denand  for  the  items  obeyed  the 
ejqxnential  distribution  and  the  holding  and  shortage  costs  were  linear, 
it  was  feasible  to  convert  tlio  integral  equation  (6.3)  into  a  hyperbolic 
partial  differential  equation  of  the  second  order  with  boundary  coiditions. 
Analytical  solution  of  tlie  boundeiry  value  problan  was  tltvJi  determined 
and  used  in  (6.1)  and  (6.2)  to  determine  Sj*,  S2*  and  C*. 


.  ^  ^  I  kw»  t  HjJ  41  1 1  ff  > 
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6.2  Extension  cind  RGccrtinciidations 

Several  interesting  directions  for  future  rcseeirch  and 
developiTont  are  suggested  as  a  result  of  this  ^research. 

The  conputationeil  asisects  of  tlie  problem  need  further  investi¬ 
gation.  The  necessary  and  sufficient  conditions,  established  in 
this  study,  for  the  existence  of  a  relative  minima  are  valid  for  any 
continuoiis  density  function  and  any  conditional  expected  holding  and 
shortage  cost  function  whicli  is  twice  differentiable.  In  tlie  case  of 
a  two-oamrxlity  problem,  where  the  demand  for  tlie  items  obeys  tlie  ej^xe- 
nential  distribution  and  tlie  holding  and  shortage  costs  are  linear, 
it  has  been  feasible  to  convert  the  integral  eqmtion  into  a  partial 
differential  equation  of  tlie  second  order  with  boundnay'  conditions. 
HoA’ever,  what  about  tlie  case  when  the  doiTiands  for  the  tivo  items  are 
given  by,  say  a  ganwa  distribution?  Is  it  tlien  feasible  to  convert  the 
integral  equation  into  a  partial  differaitial  equation?  And  if  so, 
vbat  COTputational  procedure  must  be  follwed  to  solve  for  tlie  policy 
parameters?  Another  question  also  arises,  hoi?  to  solve  for  tlie 
cptintil  policy  parameters  for  the  two-comiodity  problem,  when  the  demand 
for  the  items  are  exponentially  distributed  and  the  holding  and  shortage 
costs  are  not  linear  but,  say,  quadratic. 

As  we  have  seen  tlie  determination  of  tlie  optimum  policy  paraneters 
is  a  conplex  task  even  in  the  case  of  a  siiiplo  denond  density  function. 


IT 
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As  an  alternative  to  approximating  these  paraineters  by  numerical  iieans, 
the  approach  used  by  Roberts  [3]  could  be  followed  to  give  analytical 
formulas  for  determining  the  asynptotic  reordering  region  for  the  case 
when  the  holding  eind  shortage  costs  axe  linear.  In  order  to  derive 
the  asynptotic  e>qpressions ,  for  tlie  set  of  equations  that  is  used 
to  detemune  the  policy  parameters,  an  approximation  for  tlie  rencival 


function 


R(S,C) 


ip  (t)  ^  and  similar  extensions  of  the  type  discussed 


in  Smith  [14]  are  required.  Tlie  v;ork  of  Bickel  and  Yahav  [21]  and 
Parrel  [4]  in  investigating  the  asymptotic  beliavior  of  tlie  renewal 
function  in  two  and  higher  cinyinslo.is  must  be  noted.  Ho/zever,  the 
challenging  problem  is  to  find  an  analogous  tlieorem  to  Smith's  theoran 
involving  tlia  renewal  density  function  in  teo  and  liigher  dinensions. 

As  noted  in  Chapter  I,  under  tlie  adopted  (cr,S)  policy, 
the  sequence  of  stock  levels  at  the  beginning  of  each  period  forms  a 
discrete  Markov  process.  Greenberg's  [5]  approacli  could  be  follc^^^ed 
to  determine  the  transient  distributions  of  the  stock  levels  prior 
to  ordering.  For  tlie  stationary  distributions  Karlin  [2]  approach 
could  be  follaved.  Tlie  results  of  this  resecurch,  especially  the  proper¬ 
ties  of  the  g-convolution  operation,  are  hoped  to  be  of  great  use  in 
determining  the  transient  and  stationary  distributions  of  the  stock 
levels  prior  to  making  ordering  decision. 

Hiis  study  treated  "the  Patient  Custcansr  Cass"  where  all  vuifilled 
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demands  were  backlogned.  For  the  lost  case,  .the  analysis  can  proceed 
along  the  same  procedures  used  in  Chapters  III  and  IV.  -  , 
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p.  the  sliortage  cost  per  unit  of  cxximcxiity  i  on  hand  at  the  end 

^  of  a  period 

L(x)  the  conditional  expected  shortage  and  holding  cost  neasured 
at  the  end  of  a  period 

g{o,S)  the  stationary  total  e^^ioctcd  cost  per  period 

I^(x)  the  modified  Bessel  function  of  order  n 

starred  synbols  denote  conditions  at  optimality 

the  value  of  the  total  ej^iected  cost  per  period  exeltding 
the  variable  ordering  cost 
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NUMBER  OF  SUBDIVISIONS  N1  =  100 
OPTIMUM  VALUES  OF  C  ANO  S 
CO  =  0.y263V9993096A80  01 
ITERATION  NO.  =  4 

M(S,C)  Sl=«S2  C 

0.500069499101990  01  0 .401999664306640  01  0.826249994011600  01 
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ITERATION  NO.  =  1 

0.499293609329910  01 
ITERATION  NO.  =  I 

0.50025237931040D  01 
ITERATION  NO.  =  2 

I 

0.49061731530005D  01 
ITERATION  NO.  =  2 

0.499634547244560  01 
ITERATION  NO.  =  2 

0.499641934096120  01 
ITERATION  NO.  =  2 

0'.499768472700720  01 
ITERATION  NO.  =  2 

0.499764112277910  01 
ITERATION  NO.  =  2 

0.499701551600910  01 
ITERATION  NO.  =  2 

0.499735851909720  01 
ITERATION  NO.  =  2 

0’,499264069413920  01 
ITERATION  NO.  =  2 

0.490965975333530  01 
ITERATION  NO.  =  2 

0.490095157206750  01 
ITERATION  NO.  =  3 

0.498075712864140  01 
ITERATION  NO.  ='3 

0.490075712064140  01 
ITERATION  NC.  =  3 

0.499097963353170  01 
ITLKAVION  NO.  =  3 

0.499046233805930  01 
ITERATION  NO.  =  3 

0.50008  154  703, '0:10  01 


0.39199991226196E  01 

0.39399986267090E  01 

0.39199991226i96E  01 

0.39399986267090E  01 

0.39599901307983E  01 

0.39799976348877E  01 

0.39999971389771E  01 

0.40199966430664E  01 

0.40399966430664E  01 

0.40599956512451E  01 

0.40799951553345E  01 

0 .40999946594238E  01 

0.39199991226196E  01 

0.392g9991226196E  01 

3 .39399986267090E  01 

0.39599981307983E  01 

0.397999763488771:  01 


0.82639999389648D  01 

0.026399993896480  01 

0.826199994049970  Oi 

0.826199994049970  01 

0.826199994049970  01 

0. 82619999404997D  01 

0.826199994049970  01 

0.326199994049970  01 
0,826199994049970  01 

0.826199994049970  01 

0.826199994049970  01 

0.826199994C  '9970  01 

0.826299^93973230  01 

0.826199994049970  01 

0.826299993973230  01 

0.826299993973230  01 

0. 8262999939732:0  01 
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ITERATION  NO.  =  3 


0.500189183104950  01 
ITERATION  NQ-.  =  3 

0.500148781866140  01 
ITERATION  NO.  =  4 

t 

0.498800582182030  01 
ITERATION  NO.  =  4 

0.499763824692570  01 
ITERATION  NO.  =  4 

0.499769408400690  01 
ITERATION  NO.  =  4 

0.500003905046990  01 
ITERATION  NO.  =  4 

0.499904571616230  01 
ITERATION  NO.  =  4 

0.500069499101990  D1 
ITERATION  NO.  =  4 
0.499933712128080  01 
ITERATION  NO.  =  4 

0.499334954327350  01 
ITERATION  NO.  =  4 

0.499004532973770  01 
ITERATION  NO.  =  4 

0.490284551500160  01 


0.3999997138977] E  01 

0.40199966430664E  01 

0.39199991226196E  01 

0.39399906267090E  01 

0.39599981307983E  01 

0.39799976340877E  01 

0.39999971389771E  01 

0.40199966430664E  01 
0.40399961471550E  01 

0.40599956512451E 

0.40799951553345e  01 

0.40999946594238E  01 


0.826299993973230  01 

0.826299993973230  01 

0.02624999401160D  01 

0.026249994011600  01 

0.826249994011600  01 

0.826249094011600  01 

0.026249994011600  01 

0.826249994011600  01 
0.82624999401160D  01 

0.02624999401160D  01 

0.8262499940U600  01 

0.826249994011600  01 
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